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ABSTRACT 


A  well  graded  silt  was  subjected  to  frost  heave  tests. 

Specimens  three  inches  in  diameter  by  nine,  sixteen  and  twenty  four 
inches  long  were  frozen  from  the  top  down,  using  the  "open  system" 
of  test.  Instrumentation  was  provided  for  the  determination  of 
heave,  and  both  piezometric  pressure  and  temperature  at  intervals 
along  the  specimens.  Pressures  were  determined  by  means  of  water- 
over-mercury  manometers  connected  to  porous  tapping  points  embedded 
in  tne  speciments. 

The  variation  in  the  piezometric  pressure  along  the  flow 
path  was  found  to  be  incompatible  with  Darcy' s  law  for  the  flow  of 
water  through  a  saturated  homogeneous  material.  The  difference  was 
not  reconciled  by  the  consideration  of  variations  in  viscosity  and 
void  ratio  along  the  flow  path. 

Assuming  that  there  was  continuity  of  flow  between  the  reservoir 
and  the  frost  line,  the  permeability  of  the  silt  was  found  to  decrease 
with  increasing  negative  piezometric  pressures,  from  zero  throughout 
the  pressure  range  recorded.  This  could  not  be  explained  by  the  entry 
of  air  into  the  specimen  as  it  occurred  for  negative  pressures  that 
were  much  smaller  than  the  capillarity  of  the  material. 

The  pressure-permeability  relationships  obtained  suggested  the 
applicability  of  "unsaturated"  flow  theory  to  soil  moisture  movements 
due  to  freezing.  This  possibility  was  supported  by  similarities 
observed  in  certain  fl ow  characteristics  for  flow  produced  by  ice 
leasing,  evaporation  and  applied  negative  pore  pressures.  Furthermore, 
Gardner's  empirical  equation  for  the  pressure-permeability  curves  for 

flow  under  negative  pressure  was  fitted  to  the  frost  heave  test  data 
although  the  n-value  required  was  lower  than  any  reported  by  Gardner. 
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CHAPTER  I 


INTRODUCTION 

Certain  phenomena  related  to  the  freezing  and  thawing  of  soils 
have  been  observed  from  early  times  but  these  received  little,  if  any, 
serious  study  prior  to  the  advent  of  modern  transportation  systems.1 
Subsequently,  "frost  action"2  has  become  an  economic  problem  of  prime 
importance  to  highway,  airport  and  railroad  administrations. 

One  field  of  thought  concerning  frost  heave  stemmed  from  the 
known  fact  that  water  increases  in  volume  by  about  nine  percent  when 
it  changes  into  ice.  Between  1916  and  1930,  Taber  (l,  2,  3;  4)3 
performed  tests  which  showed  that  crystals  are  capable  of  exerting 
pressure  during  growth.  Furthermore,  he  showed  that  the  heave 
obtained  in  his  tests  was  often  greater  than  that  which  could  be 
explained  on  the  basis  of  volume  change  on  freezing,  alone.  On  the 
basis  of  his  findings,  he  suggested  that  the  difference  was  due  to 
"ice  segregation"  which  would  occur  when  there  was  migration  of 
moisture  to  the  frost  line.  This  hypothesis  was  substantiated  when 
heave  occurred  with  a  pore  liquid  which  contracted  on  freezing.  Thus, 
the  importance  of  the  availability  of  water  was  recognized  and  the 

1For  a  survey  of  past  work  in  the  field  of  frost  action  and 
related  fields,  see  "Frost  Action  in  Roads  and  Airfields  -  A  Review 
of  the  Literature  1765-1951";  A.W.  Johnson,  Highway  Research  Board, 
Special  Report  No.  1,  1952. 

2The  frost  action  terminology  used  is  in  accordance  with  the 
list  of  terms  and  definitions  approved  by  the  Highway  Research  Board 
Committee  on  Frost  Heave  and  Frost  Action  in  Soils (5).  Other 
specialized  terms  will  be  explained  as  they  occur. 

^Numbers  in  parentheses  refer  to  the  bibliographical  listings. 


i :  .. 


.  '.O  v 


• 

v  .  !  >.  . 

• 

•  •  J". 

v\  > 

«  • 

<,  J.  .  • 

• 

. 

%  ^ 

2 


possibilities  for  controlling  frost  heave  were  extended  to  the  control 
of  soil  moisture  movements. 

It  was  recognized  by  Taber,  that  the  availability  of  water  at 
the  frost  line  is  dependent  on  the  existence  of  a  source  of  supply 
within  reach  of  the  capillary  rise  of  the  material,  and  on  the  resist¬ 
ance  offered  by  the  material  to  flow  from  the  source  to  the  frost  line. 
The  first  of  these  factors  may  be  controlled  in  laboratory  tests  and 
therefore,  the  "open"  and  "closed"  system  concepts  of  Taber  must  be 
considered  when  comparing  frost  heave  test  results  from  independent 
sources. 

Based  on  both  field  and  laboratory  investigations  Taber  (4), 
Beskow(6)^,  and  Casagrande(7)  each  suggested  certain  grading  limits 
for  frost  susceptible  material.  The  Casagrande  criterion,  which  is 
commonly  used  on  this  continent,  has  very  limited  value  for  borderline 
materials  such  as  dirty  gravels.  According  to  Ducker^)-1-,  "The 
quantitative  frost  behavior  of  a  soil  is  determined  not  by  the 
percentage  of  grain  sizes  below  0.02  mm.,  but  by  the  mineralogical 
and  chemical  composition  of  the  ’binding*  material". 

Beskow(6)  performed  frost  heave  tests  on  both  natural  soils 
and  separated  fractions.  These  showed  that  the  maximum  heave  occurred 
for  grain  sizes  in  the  fine  silt  range  (using  M.I.T.  grain  size  scale). 
On  the  basis  of  investigations  into  the  effect  of  both  surcharge  and 
applied  pore  stresses  of  different  magnitudes,  Beskow  concluded  that 
"Capillary  and  load  pre s sure s  then  have  a  similar  effect  and  the 


^Reference  is  to  the  translation 
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frost  heave  can  be  stated  in  relation  to  the  total  pressure  acting 
at  the  frost  line,  which  is  the  sum  of  the  two.  The  hypothesis  used 
by  Beskow  in  determining  mathematical  expressions  relating  the  rate 
of  heave  to  the  total  pressure  at  the  frost  line,  was  that  the  pressure 
deficiency  or  negative  pore  fluid  pressure  has  a  maximum  value  equal 
to  the  capillarity1  of  the  material. 

Both  Taber (k)  and  Beskow (6)  considered  that  moisture  movements 
took  place  in  the  liquid  phase,  although  there  is  controversy  as  to 
the  significance  of  vapor  diffusion  in  these  moisture  movements.  By 
computing  the  rate  at  which  moisture  could  be  moved  by  vapor  diffusion, 
using  the  equations  given  by  Zunker,  and  comparing  with  the  rates  of 
flow  obtained  in  frost  heave  tests,  Beskow(6)  showed  that  vapor 
diffusion  could  account  for  only  about  one  thousandth  of  the  rate  of 
heave  obtained.  More  recently,  Kuzmak  and  Serda(9)>  on  the  basis  of 
studies  of  moisture  movements  in  porous  media,  concluded  that  water 
movements  due  to  suction  gradients  take  place  in  the  liquid  phase. 

The  effect  of  different  admixtures  was  investigated  by  Yurkiw(lO), 
Luck(ll),  and  Lambe(l2).  Yurkiw  showed  that  lignosol  was  effective  in 
reducing  heave,  and  that  it  reduced  the  permeability  of  the  specimens. 

He  suggested  that  the  reduction  in  permeability  was  due  to  reduced 
pore  space.  Luck  concluded  that  "The  viscosity  of  the  saturating 
fluid  is  a  contributing  factor,  but  not  a  controlling  factor,  in  the 
reduction  of  frost  heaving  of  soils".  Lambe  reported  on  a  testing 
program  in  which  a  variety  of  additives  and  soil  types  were  used. 

These  additives  were  grouped  into  the  following  classifications; 
void  pluggers  and  cements,  aggregants,  dispersants,  and  waterproof ers . 

^The  negative  pressure  required  to  draw  air  through  an  initially 
saturated  sample. 
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These  tests  showed  the  polyvalent  cations  and  the  dispersants  to 
be  the  most  promising  in  reducing  heave,  while  "dirty*'  gravels, 
sandy  clays,  and  silty  sands  were  the  most  responsive  to  treatment. 
Uniform  silts  and  lean  clays  were  the  least  responsive. 

Ruckli(l3)  and  Jumikis(l4)  gave  a  theoretical  treatment  to 
the  thermodynamical  aspects  of  the  frost  heave  problem.  Fenner (15) 
introduced  the  concepts  of  unsaturated  flow  and  of  moisture-tension 
relationships  into  the  frost  heave  problem. 

In  the  various  investigations  referred  to  above,  the  effect 
of  varying  the  length  of  the  flow  path  (rater  reservoir  to  frost  line) 
was  not  tested.  Also,  there  is  a  lack  of  data  on  the  physical 
conditions  that  exist  or  develop  in  the  pore  fluid  during  the  freezing 
process.  This  investigation  is  intended  to  provide  more  information 
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CHAPTER  II 


THEORY 

As  mentioned  previously,  frost  heave  is  due  to  ice 
segregation  which  takes  place  when  moisture  migrates  to  the 
frost  line.  Comparing  the  rate  of  moisture  movement  hy  vapor 
diffusion,  as  computed  using  equations  given  hy  Zunker,  with 
the  flow  determined  from  frost  heave  tests,  Beskow(6)  concluded 
that  movements  in  the  vapor  state  are  negligible.  More  recently, 
Kusmak  and  Sereda(9)  found  that  moisture  movements  in  porous 
material  due  to  suction  gradients  took  place  in  the  liquid  phase. 
Therefore  only  theory  of  liquid  flow  will  be  considered. 

Liquid  Flow  in  General 

It  is  characteristic  of  liquids  that  shearing  stresses 
produce  flow.  The  exi stance  of  shearing  stresses  is  dependent 
upon  differences  in  the  energy  level.  In  fluid  mechanics, 
Bernoulli's  equation  expresses  explicitly  three  energy  components 
(l6,  17)*  An  element  within  a  mass  of  liquid  will  possess  pressure 
energy  due  to  the  pressures  acting  on  that  element.  It  will  have 
elevation  or  potential  energy  by  virtue  of  its  position  with  respect 
to  any  arbitrary  datum,  and  if  the  element  is  in  motion,  it  will 
have  kinetic  or  velocity  energy.  Expressed  algebraically, 

Ev  =  V2/2g,  Ep  =  p/  ^  =  h,  and  Ee  =  z 
where  V  is  the  velocity  of  the  element,  g  is  acceleration  due  to 
gravity,  p  is  the  average  pressure  acting  on  the  element,  is  the 
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unit  weight  of  the  liquid,  h  expresses  p  in  terms  of  height  of 
liquid,  and  z  is  the  distance  from  the  element  to  an  arbitrary 
datum.  The  total  energy  is  equal  to  the  sum  of  the  three 
components. 

E  =  z  +  p/7L  +  V2/2g  =  z  +  h  +  V2/2g  (II. 1) 

The  quantities  in  equation  (il.l)  are  energies  per  unit  weight 
of  the  liquid  considered,  hence  the  units  of  length.  The  term 
h  is  meaningful  only  if  the  liquid  is  incompressible.  In  fluid 
mechanics,  liquids  are  considered  to  be  incompressible  for  steady 
flow  conditions. 

Consider  the  flow  in  a  uniform  pipe(See  Plate  l).  Provided 
that  the  various  energy  components  can  be  evaluated  at  sections 
1  and  2,  the  energy  loss  per  unit  weight  of  liquid  H,  will  be  given 
by  the  equation 

(z-l  +  h-L  +  V12/2g)  -  (z2  +  h2  +  V22/2g)  =  H 
or 

(Zl  -  z2)  +  (1^  -  h2)  +  (Vx2  -  V22)/2g  =  H  (II. 2) 

where  Vj_  and  V2  represent  average  velocities,  calculated  by  dividing 
the  volume  of  flow  per  unit  time  by  the  cross-section  area  of  the 
pipe.  The  subscripts  designate  the  sections  shown  on  Plate  1. 

A  plot  of  the  type  shown  on  Plate  2  nay  be  obtained 
experimentally  using  equation  (II.2)  and  data  from  a  pipe  flow 
apparatus.  At  the  low  end  of  the  average  velocity  scale,  the  energy 
loss  is  proportional  to  the  average  velocity,  a  condition  which 
defines  laminar  flow.  At  higher  velocities  the  flow  is  referred  to 
as  turbulent.  Here  the  energy  loss  is  approximately  proportional  to 


the  average  velocity  squared. 
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ENERGY  CONCEPTS  APPLIED  TO  FLOW  D!  A  UNIFORM  PIPE 


Average  Velocity  ■  > 


PLATE  2 


TYPES  OF  FLOW 
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Based  on  permeability  tests,  Burmister(l8)  concluded  that 
in  nature  flow  through  silt  takes  place  in  the  laminar  state. 

As  the  material  used  in  the  present  testing  program  was  a  fine 
silt,  and  because  (as  mentioned  previously)  the  optimum  grain 
size  for  producing  ice  lenses  falls  in  this  range,  only  laminar 
flow  will  be  considered. 

Viscosity  is  defined  as  "  .  .  .  that  property  of  a  fluid 
by  virtue  of  which  it  offers  resistance  to  shear  stress"  (17). 

Viscous  effects  must,  then,  contribute  to  the  energy  losses 
which  occur  in  a  flowing  liquid.  The  shearing  stress  due  to 
viscosity  is  governed  by  Newton’s  law  of  viscosity. 

=  pv/y .  (II.  3) 

where  T*  is  the  shearing  stress,  p.  is  the  viscosity,  and  v  is  the 
change  in  velocity  over  a  distance  y  normal  to  the  direction  of 
flow.  Dividing  through  by  the  unit  weight  of  the  liquid,  equation 
(II.3)  becomes 

377L  =  Hv/y/j^ .  (II.  *0 

This  is  an  expression  for  the  energy  loss  per  unit  weight  of  the 
liquid.  For  p,  y,  and  constant,  the  energy  loss  is  proportional 
to  the  change  in  velocity.  Providing  the  geometry  of  the  velocity 
distribution  curve  remains  similar,  the  average  velocity  change  will 
be  a  constant  function  of  the  maximum  velocity  change  v.  Therefore, 
the  energy  loss  per  unit  weight  of  liquid  is  proportional  to  the 
average  change  in  velocity.  If  y  is  zero  at  some  fixed  boundary 
where  the  velocity  of  flow  is  zero,  then  the  energy  loss  is 
proportional  to  the  average  velocity:  the  condition  given  previously. 
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as  that  defining  laminar  flow.  Thus,  in  laminar  flow,  the  energy 
loss  as  given  by  equation  (II.2)  is  due  to  the  viscosity  of  the 
fluid. 

Another  characteristic  of  laminar  flow  is  that  flow  lines 
(the  paths  taken  by  individual  particles  in  passing  through  a 
given  system)  do  not  cross.  It  follows  that  the  continuity 
equation  applies  to  flow  through  a  flow  tube  (the  space  bounded 
by  flow  lines).  Consequently,  the  mass  per  unit  time  passing 
through  any  section  of  a  flow  tube  is  constant  for  steady  flow 
conditions.  Thus,  knowing  the  rate  of  flow  and  the  geometry  of 
the  flow  tube,  the  average  velocity  at  any  section  can  be  determined. 

Summing  up,  it  can  be  said,  that  for  the  steady  laminar  flow 
of  an  incompressible  liquid  in  a  flow  tube,  equation  (ll.2)  may  be 
used  to  determine  the  energy  loss  due  to  the  viscosity  of  the  fluid. 

When  a  flow  tube  is  straight  and  of  uniform  size,  flow  will 
be  one -dimensional.  Furthermore,  for  the  conditions  stated  above, 
the  average  velocity  will  be  the  same  at  every  section  along  the 
flow  tube.  The  velocity  term  in  equation  (II.2)  will  be  zero,  so 
that  for  one -dimensional  flow  that  equation  becomes 

H  =  (z-l  -  z2)  +  (^  -  h2) . (II. 5) 

In  flow  equations,  energy  loss  may  be  included  in  a  term 
called  the  hydraulic  gradient.  This  is  equal  the  energy  loss  per 
unit  distance  in  the  direction  of  flow.  Thus 

i  =  H/L  . (II. 6) 

where  i  is  the  hydraulic  gradient  and  H  is  the  energy  loss  over  a 
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distance  L.  As  L  is  independent  of  factors  affecting  energy  loss, 
all  the  assumptions  governing  the  energy  loss  will  he  inherent  in 
the  hydraulic  gradient. 

Liquid  Flow  in  Pipes 

It  can  be  shown(l7)»  that  for  steady  laminar  flow  through 
uniform  circular  horizontal  pipes,  the  average  velocity  V,  is 
given  by 

v  =  h  [Vi .  (II-7) 

where  R  is  the  radius  of  the  tube,  p  is  the  viscosity  of  the  liquid, 
and  Ap  is  the  pressure  drop  over  a  distance  L.  The  derivation 
assumes  that  the  velocity  of  liquid  in  contact  with  a  boundary  surface 
is  zero.  By  substituting  p  =  h^,  from  equation  (il.l),  we  get 

V  =  .  (II. 8) 

op  L 

By  substituting  from  Equations  (II.5)  and  (II.6),  we  can  write 

V  =  (R2^8p)i .  (II. 9) 

The  terms  inside  the  brackets  are  physical  properties  of  either 

the  pipe  or  the  liquid.  This  equation  is  applicable  to  vertical  or 
inclined  pipes  as  well  as  horizontal  ones. 

Darcy1  s  Law 

In  the  field  of  soil  mechanics,  flow  equations  are  based  on 
Darcy's  law  (19,  20),  which  in  turn,  was  based  on  the  observation  of 
water  moving  vertically  downward  through  horizontal  sand  filter  beds. 

The  flow  equation  may  be  written  as  follows. 


Q  =  k  i  A  t 


(II. 10) 
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Where  Q  is  the  volume  of  flow  occurring  in  time  t  through  a  total 
cross-section  area  A,  under  the  hydraulic  gradient  i.  The  term  k 
is  the  constant  of  proportionality  in  the  equation.  Letting  the 
rate  of  flow  q  =  Q/t,  it  follows  that  q  =  k  i  A  and  q/A  =  k  i, 
where  q/A  has  the  units  of  velocity.  Designating  this  term  the 
superficial  velocity  v,  the  equation  becomes 


(II. 11) 


v  =  k  i 


where  k  is  the  Darcy  coefficient  of  permeability  and  i  is  the 
hydraulic  gradient.  In  an  apparatus  such  as  that  shown  in 
Plate  3>  'the  superficial  velocity  will  be  equal  to  the  velocity 
of  the  liquid  as  it  approaches  and  after  it  emerges  from  the  soil 
sample.  Within  the  specimen  part  of  the  total  cross-section  area 
is  occupied  by  soil  particles,  so  that  the  average  velocity  within 
the  pores  will  be  greater  than  the  superficial  velocity. 

Letting  Vv  be  the  volume  of  voids  and  V  the  total  volume 
of  the  soil  mass,  the  porosity  n  is  defined  by  (19) 


n  =  Vv/V 


(11.12) 


Assuming  that  the  distribution  of  solid  area  As,  and  void  area  Av 
on  a  plane  through  the  sample,  is  the  same  as  the  distribution  of 
solid  volume  V_  and  void  volume  V  throughout  the  specimen,  it 
follows  that  n  =  Vv/V  =  Ay/A.  Hence  Av  =  n  A.  Applying  the  concept 
of  continuity  to  the  system  gives  q  =  v  A  =  vs  Ay.  This  may  be 
rewritten  as  v  =  vs  (Ay/A).  Substituting  n  =  Ay/A  gives 


v  =  n  vs 


(n.13) 


where  v  is  the  seepage  velocity.  This  term  has  practical 
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significance  in  that  it  represents  the  rate  at  which  an  element  of 
liquid  moves  in  the  direction  of  the  flow  path.  It  may  be  considered 
an  average  velocity  based  on  the  sum  of  displacement  vectors  which 
trace  the  path  of  a  particle  through  the  specimen.  Displacement 
components  normal  to  the  direction  of  flow  are  considered  to 
balance  out  statistically.  The  significance  of  superficial  velocity 
is  that  it  represents  the  quantity  of  liquid  passing  through  a  unit 
of  cross-section  area  of  sample,  in  unit  time. 

Equation  (il.ll)  gives  Darcy's  law  in  its  simplest  and  most 
general  form.  As  stated,  the  terms  i,  and  v,  correspond  to  cause 
and  effect,  respectively.  The  physical  characteristics  of  the  soil 
and  the  permeating  liquid  are  implicit  in  k.  The  units  of  k  are 
dependent  on  the  units  used  to  express  energy  loss,  and  consequently 
on  the  units  of  hydraulic  gradient  (2l).  As  developed  in  this 
chapter,  the  hydraulic  gradient  is  unitless.  This  gives  permeability 
the  units  of  velocity. 

The  similarity  between  equation  (II.9)  for  flow  in  pipes  and 
equation  ( Il.ll)  for  flow  in  soils,  is  apparent.  Considering  the 
former  to  be  a  special  case  of  the  latter,  it  follows  that  the 
assumptions  governing  pipe  flow  will  apply  to  soil  moisture  movements. 

In  equations  (II.9)  and  (il.ll),  the  hydraulic  gradient  is 
common.  On  the  basis  that  superficial  velocity  is  proportional  to 
the  seepage  velocity,  it  can  be  accepted  that  the  average  velocity 
in  the  pipe  flow  equation  is  similar  to  the  superficial  velocity  in 

the  Darcy's  law  equation.  Finally,  the  permeability  k  corresponds 

2  i 

to  the  term  If  k  is  multiplied  by  n//^,  the  resulting 
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quantity  K  will  be  independent  of  the  permeating  liquid.  When 
the  liquid  is  water,  the  unit  weight  may  be  considered  to  be  independent 
of  pressure  and  temperature.  In  the  metric  system  it  is  equal  to  one, 
therefore,  it  may  be  ignored  numerically.  Permeability  may  be 
corrected  to  the  viscosity  at  some  standard  temperature  (j._, 
without  altering  the  permeability  units,  by  using  the  following 
equation  (19) 


k(corrected)  =  k  L  •  •  •  •  ....  .  (II.14) 

Us 

It  being  difficult  to  define  or  determine  physical  constants 
which  are  representative  of  a  soil  sample,  permeability  is 
generally  determined  directly  by  field  or  laboratory  tests.  It 
has  been  determined  by  Taylor  (19),  that  the  following  approximate 
relationship  between  permeability  and  void  ratio,  e  (volume  of 
voids  divided  by  volume  of  soil  solids)  applies  to  sands. 

-  3  * 

62  . (11.15) 


k0  =  - - 

2  1  +  e- 


1  +  e. 


It  has  been  shown  by  Lambe  (22),  that  for  kaolinite,  this  relationship 
is  not  linear.  He  explains  this  on  the  basis  of  adsorption,  a 
surface  phenomenon  associated  with  particles  of  colloidal  size. 

The  relationship  is  sufficiently  close  to  being  linear  so  that 
equation  (11.15)  may  be  used  as  an  approximation.  At  any  rate, 

Darcy’s  law  will  apply  if  the  soil  structure  remains  constant  and 
if  the  soil  is  homogeneous  and  saturated  (i.e.  the  soil  voids  are 
filled  with  water). 


Water  Flow  in  Soils 

With  ice  segregation  moisture  generally  moves  against  gravity, 
from  the  water  table  to  the  frost  line,  as  illustrated  in  Plate  4(a). 
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Elevation  Pressure  Total 

Energy  +  Energy  ~  Energy 


(c)  Energy  Distribution  for  Darcy  Type  Flow 

Upwards 


PLATE  4 


ENERGY  PROFILES  FOR  VERTICAL  SPECIMENS 
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Point  1  represents  a  free  water  surface  at  atmospheric  pressure. 

The  frost  line  at  any  given  time  is  represented  by  point  2. 

Under  hydrostatic  conditions  the  total  energy  will  not  vary 
between  1  and  2,  and  will  be  equal  to  zero  when  the  water  level 
in  the  reservoir  is  taken  as  the  datum  for  elevation  energies. 
Elevation  energy  will  increase  linearly  from  zero  at  the  reservoir 
to  a  value  L  at  the  frost  line,  where  L  is  the  vertical  distance 
between  them.  It  follows  that  the  pressure  energy  must  be  equal 
and  opposite  to  the  elevation  energy.  At  an  elevation  of  x  cm. 
above  the  free  water  surface,  there  can  be  upward  movement  of 
water  only  when  the  negative  pressure  in  the  water  at  that  point 
exceeds  -x  cm.  of  water.  Furthermore,  subject  to  the  conditions 
inherent  in  Darcy's  law,  the  change  in  total  energy  with  respect 
to  elevation  above  the  free  water  surface  should  be  independent 
of  that  elevation.  The  total  energy  profile  should  be  linear, 
starting  at  zero  at  the  free  water  surface  and  increasing  negatively 
in  the  direction  of  the  frost  line.  These  relationships  are  shown 
in  Plate  4(a)  and  (b). 

In  the  case  of  ice  segregation,  the  temperature  will  vary 
along  the  flow  path,  from  the  water  temperature  at  the  water  table 
to  freezing  temperature  at  the  frost  line.  The  viscosity,  being 
dependent  on  temperature,  will  consequently  vary  along  the  flow 
path.  It  can  be  shown,  using  equations  (il.ll)  and  (ll.l4),  that, 
other  factors  being  constant,  the  hydraulic  gradient  is  directly 
proportional  to  the  viscosity.  Furthermore,  from  previous 
consideration  of  the  hydraulic  gradient,  it  follows  that  the  total 
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energy  loss  is  proportional  to  the  viscosity.  If  the  total 
energy  loss  between  two  points  is  Hq_  when  the  average  viscosity 
between  those  points  in  then  the  total  energy  loss  H2> 
corresponding  to  any  other  average  viscosity  jig,  will  be  given  by 
the  following  equation. 

H2  =  H-.  — .  (II.16) 

This  equation  may  be  used  to  determine  the  effect  of  temperature 
variation  on  the  total  energy  profile.  This  has  been  done  in 
Plate  5*  The  striaght  line  profile  shown  in  Plate  5(a)  is  the 
type  of  profile  to  be  expected  for  Darcy  type  flow  under  uniform 
temperature  conditions.  A  linear  variation  of  temperature  with 
depth  has  been  assumed,  and  will  exist  for  steady  state  heat  flow 
through  a  homogeneous  sample.  By  taking  temperatures  from  the 
temperature  profile  and  the  corresponding  viscosities  from  the 
International  Critical  Tables,  the  viscosity  profile  was  constructed. 

The  profiles  were  then  divided  into  increments  such  that  the  viscosity 
at  the  mid-point  of  the  increment  could  be  assumed  to  be  the  average 
across  the  increment.  Using  this  viscosity,  equation  (ll.l6)  was 
solved  for  the  corresponding  total  energy  loss.  The  total  energy 
profile  corresponding  to  the  linear  temperature  variation  was  obtained 
by  taking  the  total  energy  at  the  reservoir  surface  and  adding 
successively,  the  incremental  total  energy  losses  (See  Plate  5)*  As 
elevation  energy  is  independent  of  viscosity,  changes  that  occur  in 
the  total  energy  loss  due  to  variations  in  viscosity  must  also  occur 
in  the  pressure  energy  loss,  giving  a  non-linear  pressure  energy  profile. 
If  a  total  energy  profile  is  obtained  experimentally  for  a  temperature 
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which  varies  along  the  flow  path,  it  may  he  corrected  to  a 
uniform  temperature  condition  by  applying  equation  (ll.l6). 

When  this  profile  is  linear,  Darcy's  law  is  valid. 

Pressure  energy  can  be  converted  to  pressure  by  multiplying 
by  the  unit  weight  of  water.  In  the  metric  system,  the  two 
quantities  will  have  the  same  magnitude.  In  soil  mechanics,  the 
pressure  in  soil  water  is  referred  to  as  "pore  water  pressure"  or 
"neutral  pressure".  Pressure  transmitted  from  particle  to  particle 
is  called  "effective  pressure".  The  sum  of  the  two  is  termed  "total 
pressure".  The  variation  of  these  quantities  along  the  flow  path  in 
a  homogeneous  sample  has  been  outlined  by  Rutledge(22)  for  various 
static  and  steady  flow  conditions.  Using  his  approach,  neutral, 
effective  and  total  pressure  profiles  were  constructed  for  the  case 
of  water  moving  upward  through  a  homogeneous  sample  in  accordance 
with  Darcy's  law(See  Plate  6).  It  should  be  noted  that  the  effective 
pressure  is  not  constant  throughout  the  length  of  the  sample.  Because 
effective  pressure  is  related  to  soil  density  and  void  ratio  by  way  of 
consolidation1,  there  will  generally  by  a  tendency  for  non-homogeneity 
to  develop.  Whether  effective  stresses  produce  significant  changes  in 
the  soil  structure  or  not  depends  on  the  magnitude  of  the  stress 
variation,  on  the  consolidation  characteristics  of  the  material,  and 

1  A  theoretical  treatment  of  consolidation  is  given  by  most 
standard  textbooks  on  Soil  Mechanics. 
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on  the  previous  stress  history  of  the  material.  These  factors 
have  to  he  evaluated  for  the  particular  situation  being  considered. 

Another  complication  is  inherent  in  the  existence  of  negative 
stresses  in  the  soil  voids.  Capillary  forces  will  be  present  where 
the  specimen  is  exposed  to  air.  For  a  soil  specimen  confined  by  a 
stack  of  rings,  the  perimeter  of  the  specimen  will  be  exposed  to  air. 

It  is  here  that  capillary  or  surface  forces  will  be  present  in  menisci 
in  the  soil  voids.  Negative  pressures  beyond  those  for  which  the 
menisci  are  fully  developed,  will  result  in  the  menisci  receding 
into  the  specimen,  thereby  desaturating  it. 

With  capillary  tubes  it  has  been  shown  that  the  height  of 
capillary  rise  is  inversely  proportional  to  the  diameter  of  the  tube. 
Therefore,  the  maximum  negative  stress  that  can  exist  in  the  liquid 
at  equilibrium  is  inversely  proportional  to  the  tube  diameter.  In  a 
soil  sample,  the  largest  pores  would  be  the  first  to  unsaturate, 
followed  by  others  in  order  of  decreasing  diameter.  Now  consider  a 
homogeneous  sample  in  which  the  neutral  pressure  varies  as  shown  in 
Plate  6.  If  at  some  elevation  the  stresses  are  sufficient  to  unsaturate 
the  sample,  then  the  degree  of  saturation  (volume  of  water  over  volume 
of  voids)  will  decrease  with  distance  above  that  elevation  due  to  the 
increasing  negative  pressures.  Below  that  elevation,  the  sample  will 
remain  saturated  and  Darcy’ s  law  will  apply.  For  unsaturated  flow  the 
boundaries  of  the  flow  channels  are  no  longer  defined  by  the  geometry 
of  the  soil  voids.  Part  of  the  void  space  is  occupied  by  air  and 
cannot  be  considered  available  for  water  flow.  Hence,  for  a  given 
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rate  of  flow,  the  seepage  velocity  will  increase  with  a  decreasing 
degree  of  saturation.  As  the  process  of  de saturation  is  selective 
as  to  the  size  of  pores  drained,  there  will  he  no  simple  relationship 
between  superficial  velocity  and  seepage  velocity.  For  example,  if  a 
sample  has  a  degree  of  saturation  of  fifty  percent  for  a  given  rate  of 
flow,  it  might  he  expected  that  the  seepage  velocity  would  be  doubled 
as  compared  to  the  same  rate  of  flow  in  the  saturated  condition. 

However,  other  factors  being  constant,  average  velocity  increases  as 
the  size  of  the  flow  channel  squared  (See  equation  (ll.7)«  Therefore, 
for  a  given  rate  of  flow,  the  quantity  of  flow  passing  through  channels 
which  are  saturated  when  the  overall  degree  of  saturation  is  fifty 
percent,  will  he  more  than  double  the  flow  through  the  same  channels 
when  the  degree  of  saturation  is  one  hundred  percent.  The  actual 
seepage  velocity  for  unsaturated  conditions  will  he  greater  than  the 
seepage  velocity  determined  by  applying  porosity  and  degree  of 
saturation  corrections  to  the  superficial  velocity.  The  problem  of 
determining  some  characteristic  flow  channel  dimensions  for  evaluating 
this  difference  is  analogous  to  the  problem  of  selecting  measurable 
physical  quantities  to  define  the  permeability  of  a  soil.  At  any  rate, 
the  presence  of  air  in  the  soil  acts  as  an  obstruction  to  the  water, 
and  therefore  decreases  the  permeability  of  the  soil.  For  this 
condition,  the  permeability  is  a  function  of  the  negative  stress  in 
the  pore  water.  It  is  also  a  function  of  the  moisture  content  because 
this  quantity  is  related  to  the  air  content,  to  the  degree  of  saturation, 
and  consequently,  to  the  negative  stress.  Darcy’s  law  for  the  case  of 
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unsaturated  flow  may  be  written  as  follows(l5): 

v  =  kui . (II.17) 

The  unsaturated  coefficient  of  permeability  is  designated  by  ku,  and 
the  terms  v  and  i  have  the  same  meaning  as  in  equation  (il.ll).  The 
superficial  velocity,  v  (equal  to  qj k)  is  therefore  the  volume  of 
water  passing  through  unit  area  in  unit  time.  It  follows  that  for 
steady  state  flow,  the  superficial  velocity  will  be  constant  for 
either  saturated  or  unsaturated  flow.  For  either  condition,  if  the 
seepage  velocity  and  the  hydraulic  gradient  are  known  or  can  be 
determined,  the  former  divided  by  the  latter  will  yield  a  permeability 
value.  If  this  is  repeated  for  different  hydraulic  gradients,  then, 
by  means  indicated  earlier  in  this  chapter,  it  can  be  determined  if 
saturated  or  unsaturated  conditions  prevail.  For  a  number  of  factors 
associated  with  soil  water  movement,  Plate  7  indicates  the  type  of 
variation  that  would  be  expected  for  steady  laminar  flow,  in  both  the 
saturated  and  unsaturated  condition. 
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CHAPTER  III 


APPARATUS  AND  INSTRUMENTATION 


The  frost  cabinet  for  testing  specimens  was  built  into  a  cold 
room.  The  top  of  the  specimen  was  exposed  to  the  cold  room  temperature 
while  the  bottom  was  open  to  normal  room  temperature  of  the  adjoining 
room  which  housed  the  temperature  and  pressure  measuring  apparatus. 

(See  Plate  8). 


The  Cold  Room 


The  cold  room  had  one  "walk-in"  type  of  door  and  a  smaller  one 
opening  into  the  adjoining  room.  Cooling  was  provided  by  a  Freon 
refrigeration  unit.  The  air  in  the  cold  room  was  circulated  through 
the  cooling  coils  by  means  of  fans.  Temperatures  were  controlled  by 
means  of  a  thermostat  having  a  range  of  plus  12  to  minus  35  degrees 
Fahrenheit .  The  computations  for  both  mean  and  median  temperatures 
in  the  cold  room  are  given  in  Appendix  B,  along  with  a  plot  showing 
the  variation  of  temperature  with  time.  The  following  table  summarizes 
the  results. 


TABLE  I 


TEMPERATURE  CHARACTERISTICS  OF  THE  COIN  ROOM 


Median 

Temp. 

Op 

Max. 

Temp. 

°F 

Min. 

Temp. 

Op 

Temp. 

Range 

Op 

Frequency  of 
Temp.  Cycle 
Minutes 

-7.6 

-4.0 

-9.1 

5.1 

7.8* 

+4.3 

+8.1 

+2.5 

5.6 

6.5 

+13.8 

+17.9 

+11.0 

6.9 

7.1 

*The  frequency  was  averaged  over 
at  least  three  cycles. 

Note :  Temperature  was  measured 
at  a  point  just  above  the  frost  cabinet. 
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PLATE  8 

SCHEMATIC  DIAGRAM  of 


FROST  CABINET 


27 


In  previous  experience,  no  temperature  fluctuations  were  registered 
by  thermocouples  located  one  quarter  inch  below  the  sample  surface. 
Below  this  point,  the  cold  room  temperature  fluctuations  of  the 
type  shown  in  Table  I,  will  have  no  effect. 

The  Frost  Cell 

This  apparatus  is  shown  in  Plate  9«  Detail  drawings  of  the 
various  components  are  included  in  Appendix  A. 

The  sample  container  consisted  of  a  stack  of  three  inch  diameter 
lucite  rings  which  fitted  together  as  shown.  This  over-lap  was 
designed  to  prevent  lateral  shifting  of  the  rings  during  the  placing 
of  the  sample.  A  loose  fit  was  provided  between  rings  so  that  they 
could  easily  separate  when  heaving  occurred  in  the  sample.  The  small 
drilled  hole  in  each  ring  was  for  the  insertion  of  a  thermocouple. 

The  test  specimen  was  supported  by  the  porous  bronze  plate 
in  the  water  reservoir.  The  water  level  was  maintained  above  this 
plate  so  that  water  was  made  available  to  the  sample.  Manometer 
lines  passed  through  the  bottom  of  tne  reservoir,  through  the  bronze 
plate  and  to  the  tapping  points.  The  reservoir  cap  served  to  keep 
insulation  out  of  the  reservoir  and  also  to  check  the  diffusion  of 
moisture  into  the  surrounding  insulation. 

The  water  reservoir  was  connected  to  a  water  supply  bottle 
and  a  stand  pipe,  both  in  the  equipment  room.  The  required  water 
elevation  was  marked  on  the  stand  pipe.  This  was  maintained  by 
adding  water  to  the  supply  bottle.  With  this  system  the  water  supply 
could  be  controlled  without  entering  the  cold  room  or  disturbing  the 
test  specimen. 
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The  frost  cell  (without  a  sample)  is  shown  sitting  in  its 
cradle,  on  Plate  10. 

Manometers 

The  type  of  manometer  used,  was  a  modification  of  a  simple 
water  over  mercury  U-tube.  One  branch  consisted  of  a  one  half  inch 
inside  diameter  mercury  reservoir,  while  the  other  branch  was  one 
millimeter  inside  diameter  capillary  tube  with  connections  on  top 
for  de-airing  the  apparatus  and  connecting  the  polythene  lines 
from  the  tapping  points  in  the  specimen.  The  rubber  tubing  that 
connects  the  two  branches  was  necessary  for  the  de-airing  procedure 
which  will  be  discussed  in  the  next  chapter.  The  use  of  mercury  in 
the  manometers  in  conjunction  with  the  capillary  sections  was 
designed  to  give  a  minimum  of  time  lag  between  changes  in  pressure  in 
the  sample  and  the  response  to  those  changes  by  the  manometers.  The 
length  of  the  capillary  section  was  such  that  the  water-mercury 
interface  was  contained  within  that  section  for  all  conditions  met 
during  the  testing.  Plate  11  shows  the  final  form  of  the  manometers 
used. 

Polythene  lines  were  taken  off  above  the  capillary  section  and 
connected  to  the  glass  tubing  which  passed  through  the  water  reservoir 
to  the  tapping  points.  Rubber  couplings  above  the  porous  bronze  plate 
made  the  connection  between  the  manometer  leads  and  the  tapping  points. 
The  use  of  couplings  allowed  some  adjustment  to  the  height  of 
individual  tapping  points,  facilitated  the  changing  of  points  and 
formed  a  non-rigid  connection  so  that  the  tapping  points  could  be 
lifted  if  subjected  to  frost  heaving.  Furthermore  the  de-airing  of 
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the  manometers  was  facilitated  by  doing  the  tapping  points  separately 
from  the  rest  of  the  system.  See  Plate  9  for  manometer  connections 
at  the  frost  cell. 

The  tapping  points  were  made  from  one  quarter  inch  diameter 
glass  tubing.  One  end  was  belled  out  to  one  half  inch  diameter  and 
a  tight  roll  of  fine  wire  mesh  was  fitted  into  the  throat  below 
the  bottom  of  the  bell.  The  stem  was  then  cut  off  to  leave  the 
required  length.  These  tapping  points  were  satisfactory  in 
preliminary  tests. 

Measurement  of  Heave  with  an  Extensometer 

An  Ames  dial  mounted  on  a  ring  stand  was  used  to  measure  the 
heave  in  the  test  sample.  The  shaft  of  the  dial  rested  on  a  metal 
plate  which  was  placed  across  the  top  of  the  sample  for  that  purpose. 
It  was  necessary  to  enter  the  cold  room  to  read  the  dial. 

Temperature  Measurements  by  Means  of  Thermocouples 

Temperatures  were  measured  with  a  "Brown  Electronik  Strip 
Chart  Recorder"  (Serial  534449,  Model  153X64Pl6-X-4l),  using  iron- 
constantan  thermocouples.  The  "fused  bead"  thermocouples  were 
made  of  No.  24  B  and  S  gage  wire.  The  same  wire  was  used  for  the 
leads  going  to  the  recorder. 

The  recorder  could  handle  up  to  sixteen  thermocouples.  The 
temperature  at  each  thermocouple  was  recorded  every  four  minutes. 

The  charts  had  a  time  scale  as  well  as  the  temperature  scale  with 
a  range  of  minus  50  degrees  to  plus  100  degrees  Fahrenheit. 

The  calibration  of  the  thermocouples  and  recorder  was  done 
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by  stirring  a  fifty  percent  mixture  of  ice  and  water  with  the 
thermocouples  while  the  temperature  was  being  recorded  on  the 
chart.  The  results  of  two  such  calibrations  is  given  in  Appendix  B. 
Assuming  that  the  water  was  at  32°F.,  the  recorded  temperatures 
averaged  0.5  and  1.3  degrees  too  high,  respectively.  Deviation 
from  the  average  was  small. 

The  Frost  Cabinet 

The  frost  cabinet  supported  the  frost  cell  so  that  the 
water  reservoir  protruded  down  into  an  insulated  box  which  was 
open  to  the  equipment  room  and  normal  room  temperature .  The 
top  part  of  the  frost  cabinet  was  a  shell  for  confining  loose 
vermiculite  insulation  which  was  placed  to  the  height  of  the 
test  specimen,  leaving  the  top  of  the  specimen  exposed  to  the  cold 
room  temperature.  The  leads  to  the  temperature  and  pressure 
measuring  apparatus  were  taken  from  the  sample  into  the  bottom  of 
the  frost  cabinet  and  out  into  the  equipment  room. 


:  o  ..  ■  .  -J'  .  ’  c!’>  3' '  i 


- ..  .  u  :■  .c'- 


_  .  .  .  .  . 

'  u  '.i. 


.  .  ^ .  _  . .  o  - 


O.  ■  :  ... 

* 


•■r  •  \.o 


o  J'  is  i 


.: ch,  • :  .  .  w  . . 


*  o 


cJ'  .  : 


•  o  .... 


O'  ■.  ,  [c  ,t 


*!*  •/"»«’ 


. 


.  '  u:  ■  :  .K/ni  ::..o  .  -  •  .  -_f  .  • 


.  v  ‘  -  •  . ..  / 

;  .  &i  iz 


.  v.'-  ...  •  .  .  c"  '  . 

:j  .  ..  o.  :  .  ;0  .  J  :  o'  .  •. 


:■  r.3  .o'  ■  • 


c 


..  •  . 


c.'-  . 


£■  J'J  o  J'  J  r.~  i'  ..  /  .  r;!-  i.v.  - 

■ 


CHAPTER  IV 


PROCEDURES 

Installation  and  Checking  of  Manometers 

The  glassware  was  cleaned  by  washing  with  a  saturated  solution 
of  sodium  bichromate  and  then  flushing  with  distilled  water.  When  the 
manometers  were  mounted  on  the  panel,  precautions  were  taken  to  insure 
that  they  were  plumb.  After  assembling  the  mercury  reservoir,  rubber 
section  and  capillary  tube,  mercury  was  poured  in  and  the  trapped  air 
was  worked  out.  The  level  of  the  mercury  in  the  capillary  section 
was  lower  than  that  in  the  reservoir  due  to  capillary  effects.  As 
piezometric  pressures  were  to  be  determined  using  differences  in  readings, 
capillary  effects  would  cancel  out  provided  that  they  were  constant. 

To  check  this,  the  reservoir  was  raised  by  increments  throughout  the 
heighth  of  the  capillary  section.  For  each  increment  the  difference 
in  elevation  between  the  mercury  in  the  reservoir  and  that  in  the 
capillary  section  was  determined.  The  procedure  was  repeated  at 
several  decreasing  reservoir  heights.  There  was  little  deviation 
between  values  obtained  for  the  elevation  difference  in  the  various 
capillary  sections  with  the  mercury  level  rising,  and  with  one 
exception,  the  same  was  true  when  the  mercury  level  was  falling.  This 
was  taken  to  indicate  uniformity  of  the  diameter  throughout  the  length 
of  the  capillaries.  The  average  elevation  difference  with  the  mercury 
level  rising  was  0.8  centimeters;  with  the  mercury  level  dropping, 

0.6  centimeters,  barring  the  one  exception.  The  results  of  these 
tests  are  given  in  Appendix  B. 
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Centimeter  tapes  were  used  for  determining  elevations  on 
the  manometers.  These  were  fixed  to  the  panel  with  their  zero  marks 
at  the  same  elevation.  (See  Plate  12.) 

I 

The  manometers  and  the  water  supply  system  were  then  set  up 
as  shown  in  Plate  13,  with  the  exception  that  the  tapping  points, 
rubber  couplings  and  the  porous  bronze  plate  were  removed.  The 
apparatus  was  then  ready  to  be  charged  with  water  and  de-aired. 

De-airing  the  Water  Supply  System 

After  the  water  reservoir  was  placed  in  the  frost  cabinet, 
the  water  supply  bottle  and  the  stand  pipe  were  connected  as  shown 
in  Plate  13.  The  porous  bronze  plate  was  removed  and  de-aired 
distilled  water  was  poured  into  the  reservoir  until  it  was  several 
inches  deep.  By  blocking  off  the  stand  pipe  and  applying  a  vacuum 
to  the  water  supply  bottle  a  flow  of  water  was  obtained  which  removed 
the  air  between  the  reservoir  and  the  bottle.  To  de-air  the  remainder 
the  supply  bottle  was  sboppered  and  the  vacuum  was  applied  at  the  stand 
pipe.  The  water  reservoir  was  filled  to  a  level  above  the  position 
of  the  porous  bronze  plate.  This  plate  was  de-aired  by  boiling  it  in 
water  under  a  vacuum.  It  was  then  slid  out  of  the  de-airing  dish  into 
the  water  reservoir,  where  it  was  seated. 

The  procedure  differed  for  Test  9,  where  no  effort  was  made  to 
de-air  the  porous  bronze  plate.  It  was  left  in  place  while  the  rest 
of  the  de-airing  procedure  outlined  above  was  carried  out. 

As  long  as  the  water  level  reamined  at  or  above  the  bronze 
plate,  it  was  not  necessary  to  repeat  the  de-airing  of  the  water 
supply  system  between  tests. 
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Schematic  diagram  of  Test  Apparatus 


Filling  and  De-airing  the  Manometers 


The  mercury  was  forced  up  the  capillary  tube  to  the  throat  of 
the  lateral  where  it  was  held  by  closing  the  pinch  cock  on  the  rubber 
tube  between  the  mercury  reservoir  and  the  capillary  section.  If 
the  tapping  points  and  rubber  couplings  were  in  place,  as  shown  in 
Plat  14,  these  were  removed.  This  left  the  stub  ends  from  the  manometer 
leads  sticking  up  past  the  porous  bronze  plate.  De-aired  distilled 
water  was  added  until  the  manometer  leads  were  flooded  and  the  water 
level  in  the  reservoir  was  well  above  them.  Vacuum  from  a  water 
trapped  source  was  connected  above  the  glass  cock  on  the  manometer 
board.  The  cock  was  turned  on  until  no  air  bubbles  showed  in  the 
water  flowing  toward  the  vacuum  supply.  A  visual  check  was  possible 
as  all  parts  of  the  manometer  involved  were  either  glass  or  semi¬ 
transparent  polythene.  After  closing  the  glass  cock,  the  vacuum  supply 
was  removed  and  the  pinch  cock  was  slowly  released,  allowing  the  mercury 
surface  to  settle  back  into  the  capillary  section,  drawing  water  in 
behind  it.  The  rubber  couplings  were  replaced  and  de-airing  of  the 
tapping  points  was  commenced.  This  consisted  of  boiling  them  in  water 
and  then  forcing  water  back  and  forth  through  the  roll  of  screen  in 
the  throat  until  no  air  bubbles  were  obvious.  With  the  tapping  point 
still  submerged,  both  ends  were  stoppered  and  it  was  lifted  out  and 
taken  to  the  water  reservoir  where  the  bottom  end  was  submerged  and 
the  stopper  removed.  That  end  was  pushed  into  the  submerged  rubber 
coupling  and  the  stopper  was  removed  from  the  top,  completing  the 
de-airing  of  the  manometer. 

In  Test  9  the  tapping  points  were  not  de-aired  separately. 

The  vacuum  drew  water  right  from  the  tapping  point  bell  rather  than 
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from  the  water  reservoir.  Due  to  the  resistance  to  flow  offered  by 
the  roll  of  screen  in  the  tapping  point  throat,  it  was  difficult  to 
obtain  a  strong  flow  of  water,  even  at  vacuums  approaching  one 
atmosphere. 

In  changing  from  one  test  to  the  next,  only  the  tapping  points 
had  to  be  de-aired  again.  If  frozen  in  for  a  sufficient  length  of 
time  they  would  be  pulled  right  out  of  the  rubber  couplings  by  the 
heaving  of  the  sample. 

Preparation  of  Frost  Heave  Test  Specimens 

Samples  were  wetted  down  with  de-aired  distilled  water.  With 
the  silt,  the  moisture  content  was  increased  to  something  greater  than 
the  liquid  limit  of  the  material  and  then  the  sample  was  allowed  to 
soak  for  at  least  twenty  four  hours  before  being  placed  in  the  frost 
cell.  The  gravel  samples  were  not  allowed  any  appreciable  soaking 
time. 

The  water  supply  system  and  manometers  were  charged  and  de-aired 
as  explained  previously.  The  water  level  was  lowered  until  it  was 
flush  with  the  top  of  the  porous  bronze  plate.  This  elevation  was 
marked  on  the  panel  in  the  equipment  room  by  the  water  level  stand 
pipe.  From  then  on  water  levels  were  gaged  from  that  mark. 

A  quarter  of  an  inch  of  fine  sand  was  placed  above  the  screens 
in  the  tapping  points  to  act  as  a  filter.  This  sand  was  de-aired  by 
boiling.  It  was  then  transferred  to  the  tapping  point  bells  which 
were  full  of  water.  With  the  free  water  surface  at  the  top  of  the 
bell,  the  initial  or  zero  readings  on  the  manometers  were  recorded. 

In  Test  9  the  level  of  the  mercury  in  the  reservoir  branch  and  that 
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in  the  capillary  section  were  recorded.  In  the  other  tests  the 
diameter  of  the  reservoir  was  large  compared  to  that  of  the 
capillary  tube  so  that  only  the  elevation  of  the  mercury  surface 
in  the  capillary  tube  was  recorded.  The  thermocouples  were  strung 
into  the  frost  cabinet  and  threaded  into  the  holder,  as  shown  in 
Plate  14. 

The  first  three  rings  of  the  sample  container  were  set  on 
the  porous  bronze  plate  and  the  reservoir  cap  was  placed  to  hold 
the  rings  in  position.  The  sample  was  then  placed,  using  a  rapid 
tamping  motion  with  a  one  quarter  inch  diameter  rod  to  push  the 
material  around  the  rubber  couplings  and  manometer  tapping  point 
stems.  After  the  first  three,  rings  were  added  one  at  a  time  and 
filled  in  a  similar  manner.  When  the  level  of  the  sample  came  up 
to  the  top  rim  of  a  tapping  point  the  free  water  was  taken  down  to 
the  sand  filter  and  the  specimen  material  was  tamped  inside  the  bell. 
This  procedure  was  continued  until  the  specimen  was  complete. 

The  thermocouples  were  pushed  about  half  an  inch  into  the  sample 
through  the  holes  in  the  rings  provided  for  that  purpose.  The 
extensometer  was  brought  into  contact  with  a  cover  plate  on  the  sample 
and  the  initial  reading  was  taken.  Loose  insulation  was  placed  level 
with  the  top  of  the  test  specimen.  De-aired  distilled  water  was  added 
to  the  water  supply  bottle  until  the  level  in  the  stand  pipe  indicated 
that  the  water  level  was  at  the  required  height  above  the  porous  bronze 
plate.  The  temperature  recorder,  circulation  fans  and  cooling  unit  were 
started  and  the  starting  time  was  recorded.  The  frequency  of  taking 
readings  varied.  Generally,  all  readings  were  taken  once  an  hour  during 
the  day  and  not  at  all  during  the  night,  except  for  temperatures  which 
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were  recorded  continuously. 

Freezing  of  Samples 

Samples  were  frozen  until  the  frost  line  was  in  the  vicinity 
of  the  highest  manometer  tapping  point,  at  which  time  they  were 
allowed  to  thaw.  This  procedure  was  followed  in  tests  9,  10,  11, 

12  and  13,  with  the  exception  of  test  11  which  thawed  considerably 
during  test  due  to  a  power  failure.  Tests  12A-12F  and  13A-13C  are 
refreeze  tests  on  the  samples  from  tests  12  and  13,  respectively. 
Samples  were  completely  thawed  and  the  sample  container  rings  were 
reseated  before  each  refreeze  test.  In  test  13,  insulation  was 
removed  from  around  the  upper  end  of  the  specimen  to  force  the  frost 
line  down;  in  other  tests  either  the  temperature  in  the  cold  room 
was  lowered  or  sufficient  time  was  allowed  for  the  frost  line  to 
penetrate  under  constant  room  temperature.  Plate  15  is  a  photograph 
of  the  specimen  from  test  13  with  the  insulation  taken  away  while  it 


was  still  frozen. 
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FROZEN  SPECIMEN 


CHAPTER  V 


RESULTS  AND  DISCUSSION 


Sample  Materials 


The  materials  used  in  the  frost  heave  tests  were  a  silt  and 
two  gravels.  The  data  and  results  of  classification  tests  on  these 
materials  are  included  in  Appendix  B.  The  grain 'size  curves  are  on 
Plate  l6,  and  the  pressure-void  ratio  curve  for  the  silt  is  on 
Plate  17. 

The  following  is  a  summary  of  the  results  of  these  tests  on 
the  silt  material. 


A.  Atterberg  Limits 


1.  Liquid  limit 

• 

25.2# 

2.  Plastic  limit 

• 

21.5# 

3.  Plasticity  index 

• 

3»7# 

B.  Hydrometer  Analysis 

1.  Sand  sizes1- 

• 

30$ 

2.  Silt  sizes 

• 

62 £ 

3.  Clay  sizes 

• 

8 i 

4.  Dpo  size  (mm.)  . 

• 

0.003 

5.  Uniformity  coefficient  . 

• 

10 

C.  Consolidation  Test 

1.  Compressive  index  .  . 

• 

0.066 

2.  Permeability  (cm. /sec.) 2 

• 

10“° 

D.  Capillarity  (cm.  of  water)  . 

• 

605 

The  grain  size  curve  indicates  that  the  material  is  a  well  graded 
silt.  The  Atterberg  limits  indicate  slight  plasticity.  On  a 


plasticity  chart,  the  material  classified  as  an  "ML"  which  includes 
inorganic  silts,  rock  flour,  etc.3.  Table  XV,  Appendix  B,  is  a 


1  Grain  sizes  are  according  to  the  Massachusetts  Institute 
of  TechnolOQ'-  grain  size  scale. 

2 

Approximate  order  of  values  computed  from  consolidation  test 
results  for  a  remolded  specimen. 

3  Arthur  Casagrande,  Classification  and  Identification  of  Soils, 
American  Society  of  Civil  Engineers,  Trans.,  Vol.  113,  1948,  P»  910 . 
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summary  of  permeability  computations.  These  were  computed  for  four 
increments  of  pressure  between  O.O69  and  l.l8  kg. /cm.  The  average 

void  ratio  in  individual  increments  varied  from  0.60  to  0.5*4  and  the 
permeability  values  ranged  from  2.4  x  10"^  to  2.0  x  10"  ^  cm. /sec.,  for 
the  first  loading  cycle.  The  corresponding  void  ratio  variation,  for 
the  first  reloading  cycle  was  0.54  to  0.53;  for  the  second  it  was  the 
same,  to  two  significant  figures.  For  the  first  loading  cycle,  the 
permeability  at  the  higher  void  ratio  was  1.2  times  that  at  the  lower 
void  ratio.  This  factor,  for  the  reloading  cycles  (determined  using 
equation  (11.15)),  is  I.05.  These  results  are  based  on  a  consolidation 
test  on  a  specimen  consisting  of  the  sixt  material  used  in  the  testing 
program.  As  the  consolidation  specimen  was  molded  at  a  moisture  content 
approximately  the  same  as  that  at  which  the  frost  heave,  specimens  were 
formed,  it  is  proposed  that  as  an  approximation  the  consolidation 
results  can  be  applied  to  the  frost  heave  specimens,  within  the  range 
of  pressures  indicated. 

The  following  is  a  summary  of  the  sieve  analysis  results  for 
the  gravels. 


Calgary  Gravel 

1.  Gravel  sizes 
Sand  sizes 
Finer  than  sand 


2. 

3. 

4. 

5. 


J10 

Uniformity  coefficient 


B.  Calgary  Gravel  (with  excess  of  material  passing  No. 4  sieve) 


1.  Gravel  sizes 

2.  Sand  sizes 

3.  Finer  than  sand 

4 •  .... 

5.  Uniformity  coefficient 


72^ 

27# 

0.5  mm 

50 


62# 
36. 5# 
1.5# 
0.4  mm 
43 


Both  of  these  gravels  classify  as  "gravels  with  fines"  (GF)  according 
to  the  Airfield  Classification  System,  and  as  "coarse  material  with 
fines"  (A-2),  using  the  Public  Roads  Classification.  The  Casagrande 
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criterion  for  frost  susceptibility  in  a  well  graded  material  is  that 
the  fraction  finer  than  0.02  mm.  shall  be  in  excess  of  three  percent  (7). 
On  this  basis  neither  of  these  gravels  should  be  frost  susceptible. 

Specimen  Data 

In  this  investigation  five  specimens  were  instrumented  and 
subjected  to  freezing.  Three  were  silt  specimens  having  different 
lengths  and  two  consisted  of  gravel  overlying  the  same  silt.  Specimen 
data  and  details  of  the  instrumentation  are  given  in  Plates  l8A  to 
l8D.  The  designation  used  for  thermocouples  and  manometers  corresponds 
with  that  used  on  the  original  data  sheets.  The  lengths  giving  the 
locations  of  thermocouples  and  manometer  tapping  points  are  accurate 
to  +  0.05  inches. 

Frost  Heave  Test  Data  and  Initial  Computations 

The  quantities  measured  during  the  frost  heave  tests  consisted 
of  the  elevation  of  the  top  surface  of  the  specimen  relative  to  an 
arbitrary  datum,  the  elevation  of  the  water-mercury  interface  in  each 
manometer,  the  time,  and  the  temperature  at  various  intervals  through¬ 
out  the  length  of  the  specimen. 

r 

Top  surface  elevation  was  determined  by  means  of  an  extensometer 
as  shown  in  Plate  15.  This  was  set  and  read  before  temperature  was 
dropped.  Subsequent  readings  indicated  that  the  surface  dropped 
slightly  before  starting  a  prolonged  rise.  The  initial  drop  is  attributed 
to  the  difference  in  contraction  between  the  sample  container  and  the 
stand  supporting  the  extensometer.  The  minimum  reading  was  taken  as 
the  initial  value.  The  difference  between  this  and  later  readings  was 
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Distance  Below 


Initial  Surface  Elevation  Initial  Surface 


PLATE  ISA 


DETAIL  OF  INSTRUMENTATION 
FOR  TEST  9 


50 


#11  Thermocouple  Distance  Below 


Initial  Surface  Elevation  i  Initial  Surface 


DETAILS  OF  INSTRUMENTATION 
FOR  TESTS  10  AND  11 


r>I 


#11  Thermocouple 

Sample  silt  -  Soil  (2) 


#10  Distance  Below 


PLATE  18C 


DETAILS  OF  INSTRUMENTATION 
FOR  TEST  12-12F 


DETAIIS  OF  INSTRUMENTATION 
FOR  TEST  13-13C 
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taken  as  the  total  heave. 

Water-mercury  interface  elevations  were  read  directly  from 
centimeter  tapes  on  the  manometer  hoard.  The  elevation  corresponding 
to  atmospheric  pressure  was  determined  as  outlined  in  Chapter  IV. 

The  difference  between  the  initial  and  subsequent  elevations  was  used 
to  determine  the  pressure  at  each  tapping  point.  The  method  used 
for  determining  pressures  is  shown  on  Plate  19.  These  quantities 
were  expressed  in  "cm.  of  water"  but  would  be  numerically  the  same, 
if  expressed  in  the  units  "gra./cm.^".  The  first  units  correspond  to 
"pressure  energy"  or  "pressure  head"  in  fluid  mechanics  terminology. 

With  the  second  units,  the  quantities  correspond  to  "pressure"  in 
fluid  mechanics,  and  "neutral  pressure"  in  soil  mechanics.  In  any 
case,  the  pressure  in  the  fluid  at  a  particular  tapping  point, 
relative  to  atmospheric  pressure,  is  specified.  These  quantities 
will  be  called  "piezometric  pressures". 

In  determining  the  elapsed  time  from  the  start  of  test,  zero 

\ 

time  was  taken  at  the  instant  the  cooling  unit  was  started. 

By  the  process  of  interpolation,  the  position  of  the  32°F. 
isotherm  was  determined.  It  was  assumed  that  the  distance  between 
thermocouples  was  the  same  as  measured  before  the  start  of  the  test 
and  that  the  temperature  gradient  was  constant  between  the  thermocouples 
considered.  The  frost  line  was  considered  to  exist  at  the  location  of 
the  32°F.  isotherm. 

All  the  original  data  is  on  file  with  the  Department  of  Civil 
Engineering,  University  of  Alberta,  along  with  quantities  computed 


directly  from  that  data 
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A  General  Consideration  of  Frost  Heave  Test  Results 

Heave,  depth  of  frost  line  and  piezometric  pressures  have  been 
plotted  as  functions  of  time,  on  Plates  20A  to  20E.  The  plots  present 
an  overall  picture  of  the  behavior  of  these  quantities  throughout 
each  test. 

During  the  course  of  freezing,  negative  pressures  developed 
in  the  pore  water.  In  certain  instances  the  pressure  values  were 
subject  to  fluctuations.  In  Tests  9  and  11  where  the  cold  room 
temperature  was  held  constant  for  extended  periods  of  time,  the 
fluctuations  were  most  pronounced.  Test  9  from  320  and  780  hours 
illustrates  the  course  of  events  following  a  drop  in  the  cold  room 
temperature.  Three  stages  can  be  identified.  In  the  first,  heave, 
depth  of  frost  line  and  piezometric  pressure  increase  relatively 
rapidly.  Second,  values  seem  to  oscillate  about  a  steadily  increasing 
average.  The  third  stage  is  characterized  by  unpredictable  fluctuations 
terminating  with  a  near  static  frost  line,  much  reduced  piezometric 
pressures  and  very  little  heave.  The  conditions  in  this  stage  seem  to 
be  approaching  equilibrium  for  the  temperatures  imposed.  Under  conditions 
where  the  average  applied  freezing  temperature  decreased  slowly  and 
steadily,  it  would  be  expected  that  stage  one  and  three  would  not  exist. 
Therefore  stage  two  conditions  would  be  expected  to  exist  in  the  field 
where  the  frost  line  is  not  affected  by  daily  air  temperature  fluctuation. 
VHaen,  over  a  period  of  time,  the  fluctuating  quantities  can  be 
represented  by  increasing  straight  lines,  the  conditions  would  appear 
to  be  steady,  with  the  average  heave,  depth  of  frost  line  and 
piezometric  pressures  all  increasing. 
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An  increase  in  the  depth  of  the  frost  line  will  result  in  a 
shorter  path  from  the  reservoir  and  a  larger  surcharge  at  the  frost 
line.  It  has  been  shown  by  Beskow  (6)  that  the  rate  of  heave 
varies  inversely  as  the  effective  pressure  at  the  frost  line;  the 
effective  pressure  being  equal  to  the  negative  pore  pressure  plus 
the  surcharge.  Therefore,  the  dropping  frost  line  results  in 
factors  whose  effects  on  the  rate  of  heave  tend  to  be  compensating. 
Increasing  negative  stresses  within  the  specimen  will  result  in 
increasing  average  hydraulic  gradients  (the  piezometric  pressure 
at  the  water  reservoir  being  constant)  and  increasing  effective 
pressures  at  the  frost  line.  These  factors  will  tend  to  produce 
opposite  effects  in  the  rate  of  heave,  assuming  that  the  flow  is 
governed  by  Darcy's  law.  For  un saturated  conditions  and  steady  flow, 
the  effect  on  the  rate  of  heave  of  increasing  the  hydraulic  gradient 
will  be  opposed  by  a  corresponding  decrease  in  the  permeability.  The 
indications  are  that  while  there  are  variations  in  a  number  of  factors 
during  the  freezing  process,  the  changes  in  individual  factors  have 
opposite  effects  on  the  rate  of  heave.  For  a  constant  rate  of  heave 
the  net  effect  is  zero.  This  situation  is  approached  on  an  average 
value  basis,  by  the  conditions  in  stage  two  (described  previously). 

According  to  the  "fluctuating  frost  line"  theory  of  Benkleman 
and  Ohmstead  (23),  an  oscillating  frost  line  is  a  prerequisite  to  ice 
segregation.  Results  indicate  that  fluctuation  of  the  ftost  line  is 
but  one  of  the  consequences  of  ice  segregation.  The  Benkleman  and 
Ohmstead  theory  assumed  that  the  frost  line  would  rise  and  fall  with 
the  surface  temperature.  While  this  may  well  be  true,  test  results 


show  that  the  frost  line  can  fluctuate  while  the  end  temperatures 
are  held  constant.  It  appears  that  the  same  situation  would  exist 
for  slowly  falling  surface  temperatures. 

According  to  the  "rhythmic  ice  banding"  theory  of  Martin  (24), 

".  .  .  the  problem  involves  the  interrelationship  of  (a)  the  phase 
change  (liquid  to  solid),  (b)  the  mass  transport  of  liquid  to  the 
freezing  front,  and  (c)  the  condition  of  a  general  unsteady  heat  flow 
found  in  a  freezing  soil."  He  goes  on  to  explain  the  cyclic  nature 
of  ice  segregation.  "Nucleation  of  ice"  represents  the  start  of  a 
new  lens.  This  takes  place  below  the  existing  frost  line.  Nucleation 
is  followed  by  a  period  of  rapid  crystal  growth,  after  which  the  growth 
terminates.  Subsequently  the  cycle  repeats.  The  physical  conditions 
during  this  cycle,  were  explained  by  Martin  as  follows.  Nucleation 
takes  place  at  some  temperature  lower  than  the  normal  freezing 
temperature  of  the  liquid.  The  rapid  growth  following  nucleation 
involves  the  freezing  of  the  "free"  water  which  is  immediately 
available.  This  is  accompanied  by  a  rapid  release  of  heat  at  the 
growing  crystal  and  by  the  depletion  of  the  moisture  in  the  vicinity. 

The  latter  factor  results  in  increasing  negative  stresses  in  the  pore 
fluid  and  lowers  the  freezing  point  temperature.  The  net  effect  of 
the  various  factors  is  to  reduce  the  rate  of  crystal  growth.  Eventually 
the  temperatures  below  become  depressed  to  a  point  where  nucleation 
occurs  and  the  whole  cycle  repeats.  It  would  appear  that  the 
fluctuations  in  heave,  depth  of  frost  line  and  piezometric  pressure 
values  observed  in  the  frost  heave  tests  can  be  explained  in  terms 
of  the  "rhythmic  ice  banding"  theory. 


On  the  basis  of  laboratory  tests  using  an  apparatus  in  which 
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the  formation  of  ice  lenses  could  be  observed,  Daxelhofer  (25) 
concluded  that  "The  formation  of  successive  ice  lenses  appears  to 
depend  on  thermal  variations  which  bring  about  thawing  at  the 
boundary  of  the  frozen  zone  ..."  This  hypothesis  differs  from 
the  theory  of  Benkleman  and  Ohmstead  in  that  the  thermal  variations 
are  not  attributed  to  variations  in  the  applied  temperatures. 

Martin  does  not  include  thawing  as  part  of  the  cycle  in  the 
"rhythmic  ice  banding"  theory,  although  thermal  variations  are 
recognized  as  "the  condition  of  general  unsteady  heat  flow".  The 
most  obvious  thawing  that  occurred  during  the  present  investigation, 
took  place  when  the  average  rate  of  heave  and  the  average  rate  of 
penetration  of  the  frost  line  were  approximately  zero.  Thus, 
thawing  and  refreezing  at  the  frost  boundary  does  not  necessarily 
result  in  the  formation  of  an  ice  lens.  In  other  cases  where  the 
piezometric  pressures  were  fluctuating  the  rate  of  penetration  of 
the  frost  line  appeared  to  approach  zero  at  the  peak  of  the 
osciHations  but  in  most  cases  there  was  no  indication  of  thawing. 
However,  the  position  of  the  frost  line  was  determined  by  an 
interpolation  in  which  it  was  assumed  that  the  distance  between 
thermocouples  was  the  same  as  at  the  time  they  were  placed  and 
that  the  temperature  gradient  was  constant  over  the  interval 
interpolated.  It  is  possible  that  these  assumptions  could  have 
introduced  sufficient  error  to  mask  minor  thawing  in  the  specimen. 
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LEGEND  FOR  PLATES  20A  TO  20E 


Tapping  Point 
Elevation  Cm. 

26.6 

22.6 

Test  9  17.7 

12.7 

33.1 

27.3 

Test  10  22.5 

17.9 

13.3 

33.1 

27.3 

Test  11  22.5 

17.9 

13.3 

12.5 

10.0 

Test  12  7.5 

5.1 

2.6 

12.5 

10.0 

Test  13  7.5 

5-1 

2.6 


Designation 
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Analysis  of  the  Pressure  Data 

The  total  head  (i.e.  total  energy  per  unit  weight  of  liquid) 
values  computed  using  equation  (il.l)  are  listed  with  the 
corresponding  piezometric  pressure  values,  in  Tables  III(a)  co  III(e). 
Plates  21A  to  21 J  show  total  head  profiles  at  selected  times  in  the 
various  tests.  In  the  early  stages  of  test  the  total  heads  were 
positive,  indicating  that  there  is  movement  of  water  downward. 

These  values  decreased  from  the  top  downward  through  the  sample  until 
total  heads  were  all  negative,  indicating  the  movement  of  water  from 
the  water  reservoir  to  the  frost  line. 

The  negative  total  head  profiles  display  curvature  of  varying 
degrees.  Some  curvature  is  to  be  expected  due  to  the  fact  that  the 
viscosity  increases  along  the  flow  path.  Viscosity  corrections  were 
made,  using  equation  (ll.l6)  and  the  method  illustrated  in  Plate  5 
(Chapter  II ). 

Tables  IV(a)  to  IV(e)  include  the  temperature  data  corresponding 
to  each  of  the  total  head  profiles  in  Plates  21A  to  21 J.  The 
temperature  profiles  (Plates  22A  to  22H)  were  used  in  conjunction  with 
viscosity- temperature  data  from  the  International  Critical  Tables 
(Plate  23)  to  determine  the  viscosity  of  the  water  at  each  manometer 
tapping  point.  One  total  head  profile  from  each  of  Test  10,  11,  12B, 
and  13B  was  corrected  to  the  viscosity  at  20°C.  The  computations  are 
listed  in  Table  V.  Total  head  profiles  in  Plate  2h  show  that  for 
constant  viscosity  conditions  the  distribution  of  total  head  along  the 
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flow  path  is  non-linear.  It  remains  to  investigate  the  homogeneity 
of  the  specimens. 

In  order  to  evaluate  consolidation  effects  in  the  test 
specimens  it  is  necessary  to  determine  the  effective  stresses. 

For  the  frost  heave  test  specimens  the  effective  stress  will  he 
equal  to  the  surcharge  due  to  overlying  material  and  apparatus 
(a  positive  quantity)  minus  the  neutral  pressure(a  negative  quantity). 
The  surcharge  pressure  per  unit  length  of  specimen  was  calculated  in 
Table  VI,  using  estimated  values  for  the  void  ratio  and  the  degree  of 
saturation.  The  maximum  error  due  to  these  assumptions  was  determined 
to  be  +5 $.  The  surcharge  pressure  at  each  tapping  point  is  listed  in 
Table  VII.  Having  piezometric  pressures  and  surcharge  pressures,  the 
effective  pressures  were  computed  (See  Tables  III(A)  to  111(e)).  The 
corresponding  void  ratios  were  taken  from  the  consolidation  curve, 
Plate  17. 

An  approximate  relationship  between  permeability  and  void  ratio 
is  given  by  equation  (II.15),  k^_  :  kg  =  +  el):e2^/ (l  +  ^2)' 

For  Darcy  type  flow  the  permeability  may  be  expressed  k  =  v/i  =  Lv/H 
(equation  (il.ll)).  For  a  particular  choice  of  time  and  distance 
along  the  flow  path,  v  and  L  are  constants,  so  that  k  is  inversely 
proportional  to  the  head  loss  H.  Substituting  into  equation  (II.15), 
we  get 


(v.18) 


H 


1 


Using  this  equation  and  the  viscosity  corrected  total  head  values, 
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total  head  values  corrected  for  both  viscosity  and  void  ratio  were 
obtained(See  Table  VIII).  Total  head  profiles  with  the  various 
corrections,  are  shown  in  Plate  24.  With  both  viscosity  and  void 
ratio  changes  accounted  for,  the  curvature  in  the  profiles  remains; 
a  situation  which  is  incompatible  with  steady  laminar  flow  in  a 
homogeneous  saturated  specimen. 

Since  the  specimens  were  formed  by  rodding  successive  layers, 
it  is  to  be  expected  that  some  air  might  be  trapped  in  the  process. 

On  the  basis  of  the  consolidation  specimen,  the  degree  of  saturation 
of  specimens  at  the  time  of  molding  should  have  been  in  the  order  of 
9 6$.  This  corresponds  to  0.035*+  cm.  3  of  air  to  1  cm.  3  of  water.  If 
the  air  was  present  as  isolated  bubbles  through  out  the  specimen,  the 
air  pressure  would  be  equal!  to  the  pressure  in  the  surrounding  fluid. 

By  Henry's  law1,  the  solubility  of  a  gas  in  a  liquid  is  directly 
proportional  to  the  pressure.  As  the  neutral  pressures  decreased 
during  the  course  of  the  frost  heave  tests,  according  to  Henry's  law, 
air  should  be  released  from  the  pore  water,  provided  that  it  was 
saturated  at  the  beginning.  The  maximum  negative  stress  of  6k2  cm.  of 
water  was  recorded  in  Test  13B  at  110.2  hours  (Table  111(e)).  The 
corresponding  absolute  pressure  would  be  1033  -  642  =  391  cm.  of  water. 
Assuming  that  the  trapped  air  was  at  atmospheric  pressure  (1033  cm.  of 
water,  absolute)  and  that  the  pore  water  was  air  saturated  initially, 
a  pressure  drop  to  391  cm.  of  water  should  result  in  the  dissolved  air 
content  dropping  to  391/1033  c»r  O.38  of  its  initial  value.  At  a 
pressure  of  one  atmosphere  and  at  20°C. ,  1  cm. 3  of  water  can  retain 

"''See  any  standard  textbook  on  the  Fundamentals  of  Physical  Chemistry. 
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approximately  O.OI87  cm.^  of  dissolved  air  (26).  At  a  pressure 
of  391  cm.  of  water,  O.OO707  cm. 3  could  be  retained.  Thus  0.0116  cm.  3 
would  be  released,  increasing  the  undissolved  air  content  to 
0.0470  cm.  /cm.  .  However,  the  effect  of  temperature  variation  has 
yet  to  be  included.  Samples  were  molded  at  room  temperature,  which 
for  the  present  purposes  will  be  assumed  to  be  20°C.  During  test 
temperatures  were  dropped  throughout  the  system.  For  the  test,  time 
and  elevation  indicated  above,  the  temperature  was  approximately  0°C. 
Between  20  and  0°C.  the  saturated  air  capacity  of  water  changes  frpm 
O.OI87  to  0.0292  cm.  3/ cm.  3,  so  that  0.0105  cm.3/cm.3  could  be  absorbed, 
leaving  the  degree  of  saturation  and  the  quantity  of  undissolved  air 
approximately  the  same  as  when  the  sample  was  molded.  Under  conditions 
imposed  during  test  changes,  the  solubility  of  air  in  water  would 
appear  to  be  small  when  the  trapped  air  is  isolated  from  the  atmosphere. 
Where  a  specimen  becomes  unsaturated  due  to  negative  pore  stresses 
exceeding  the  capillarity  of  the  material,  the  air  space  in  the  voids 
will  be  vented  to  the  atmosphere.  Air  absorbed  by  the  water  due  to 
temperature  decrease  will  be  replaced  from  the  atmosphere  and  the  degree 
of  saturation  of  the  specimen  will  be  unaltered.  Hence  the  solubility 
factor  can  be  ignored. 

When  negative  stresses  in  the  frost  heave  specimens  exceed  the 
capillarity  of  the  specimen  material,  it  is  to  be  expected  that  there 
will  be  curvature  in  the  total  head  profiles  (See  Chapter  II ). 

Variations  in  hydraulic  gradient  and  permeability  due  to  increasing 
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negative  stresses  were  predicted  as  indicated  on  Plate  7.  These 
quantities,  computed  from  experimental  data,  are  listed  in  Table  IX 
and  plotted  on  Plates  25  and  26.  The  experimental  hydraulic  gradient 
profiles  are  similar  to  those  predicted  on  the  basis  of  capillary 
desaturation.  A  similar  comparison  of  the  permeability  profiles 
indicates  that  the  permeability  decreases  in  the  direction  of  flow 
as  predicted  but  with  reversed  curvature.  The  curvature  of  the 
permeability  profile  is  dependent  on  the  hydraulic  gradient  profile. 

The  permeability  times  the  hydraulic  gradient  must  equal  the 
superficial  velocity,  which  is  a  constant  along  the  flow  path,  for 
steady  flow  conditions.  It  can  be  shown  that  the  curvature  of  the 
permeability  profile  can  be  reversed  by  increasing  or  decreasing 
the  curvature  of  the  hydraulic  gradient  profile.  Thus,  the 
experimental  permeability  profiles  are  similar  to  those  predicted. 

It  was  shown  previously,  that  the  pressures  in  the  pore  water 
increases  negatively  from  zero  at  the  reservoir  to  a  maximum  at 
the  frost  line.  If  the  capillarity  of  the  material  falls  within 
this  range  of  pressures,  it  would  be  expected  that  the  test  specimen 
would  contain  one  zone  having  a  constant  degree  of  saturation  near 
100$,  and  a  second  in  which  the  degree  of  saturation  decreases  with 
increasing  negative  pressure.  If  water  moves  through  the  saturated 
zone  subject  to  Darcy's  law,  it  should  be  registered  on  the  permeability 
and  hydraulic  gradient  nrofiles  as  vertical  lines,  indicating  that 

t 

those  quantities  are  constants  within  that  zone.  On  this  basis,  no 
saturated  zone  can  be  identified  (Plates  25  and  26). 


, 

. 

- 

.  J  X-  '  :  .'.I..'- 

-  .  . 

■  .  (  u  .  -  •.  '■  ■'  .  • .  ,  . 

■ 

. 

:  :p 


. 

[i  •  -  •  •  ' 


. 


. 


. 

. 


711 


One  capillarity  determination  was  carried  out  using 
apparatus  similar  to  Beskov's  simplified  capillarimeter  (6), 
with  the  exception  that  a  water  over  mercury  column  was  used  to 
apply  negative  stress  to  the  specimen.  The  capillarity  value 
obtained  was  605  cm.  of  water.  This  falls  within  the  range  of 
400  -  1000  cm.  given  by  Beskow  for  fine  silt  in  the  loose  state. 

It  seems  reasonable,  considering  that  the  specimen  material  was 
a  well  graded  silt  having  a  D^q  value  in  the  fine  silt  range.  If 
the  determined  value  of  capillarity  is  correct,  then  the  only  place 
that  desaturation  should  have  occurred  is  in  the  vicinity  of  the 
top  manometer  for  Test  13B,  at  107.8  to  110.2  hours.  This  indicates 
that,  while  the  profiles  discussed  above  displayed  characteristics 
that  could  be  explained  for  negative  stresses  greater  than  the 
capillarity  of  the  material,  in  fact,  stresses  were  less  than  the 
capillarity.  For  water  movement  during  ice  segregation,  it  appears 
that  the  hydro  dynamical  characteristics  of  flow  are  not  explained 
by  the  theory  presented  in  Chapter  II.  Further  evidence  to  support 
this  finding  is  to  be  found  on  Plate  27  which  shows  that  the 
permeability  decreased  with  increasing  negative  pressures  throughout 
the  measured  pressure  range  and  that  corrections  accounting  for 
changes  in  temperature  and  void  ratio  did  not  alter  the  situation. 
Pressure-permeability  curves  are  shown  for  various  times,  on  Plate  28 
which  was  plotted  from  the  data  in  Table  X.  Although  the  curves  tended 
to  shift  to  the  left  with  time,  they  retained  a  similarity  of  shape. 
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The  most  of  the  shifting  tools:  place  during  the  first  thirty- 
four  hours  of  the  test.  Plate  20E  shows  that  this  was  the 
period  in  which  the  initial  drop  in  the  cold  room  temperature 
produced  relatively  rapid  increases  in  the  heave,  depth  of 
frost  line  and  the  piezometric  pressure.  The  thirty  four  hour 
pressure -permeability  curve  appears  to  be  coincident  with  the 
subsequent  curves.  This  suggests  that  apart  from  effects  due 
to  abrupt  changes  in  the  applied  freezing  temperature,  the 
pressure -permeability  relationship  is  defined  by  a  unique  curve. 

The  sensitivity  of  the  permeability  to  pressure  changes 
increased  with  decreasing  pressure  values.  Extreme  sensitivity 
was  exhibited  as  the  negative  pressures  approached  zero.  This 
appears  to  dismiss  any  explanation  of  the  pres sure -permeability 
behavior  based  on  the  assumption  that  the  capillarity,  as  determined, 
was  not  a  reliable  measure  of  the  minimum  stress  required  to  induce 
air  entry  in  the  specimen.  Some  explanation  other  than  that  outlined 
in  the  theory  must  be  sought. 

Flow  Under  Negative  Pressure  Conditions 

It  was  noted  by  Beskow  (6)  that  the  movement  of  water  to  a 
frost  line  was  similar  to  the  movement  of  water  to  a  surface  where 
evaporation  was  taking  place.  In  either  case,  water  in  the  liquid 
phase  was  removed  from  the  soil  by  the  process  and  the  resulting 
moisture  deficiency  induced  the  subsequent  flow. 

Plate  29  shows  the  hydraulic  gradient  profiles  at  various 
times  during  Test  10.  Richards,  Gardner  and  Ogata  (27)  have 
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published,  a  similar  plot  for  the  case  of  evaporation  following 
the  prolonged  irrigation  of  a  field  plot.  This  is  shown  on 
Plate  30.  The  similarity  between  the  hydraulic  gradient  profiles 
for  the  two  cases,  indicates  a  corresponding  similarity  in 
hydrodynamics  involved  in  each.  For  example,  it  follows  that 
piezometric  pressure  and  total  head,  in  the  evaporation  case, 
will  vary  in  a  similar  manner  to  that  determined  from  the  frost 
heave  tests.  The  data  of  Richards  and  Weeks  (28)  are  from  tests 
in  which  water  was  drawn  from  a  soil  column  by  applying  tension 
to  the  pore  water  at  one  end,  after  the  soil  and  moisture  had  come 
to  equilibrium  under  a  lesser  tension.  On  Plate  JL  is  a  comparison 
between  the  data  from  the  frost  heave  tests  and  the  equivalent  values 
for  the  case  of  applied  tension.  The  variation  in  the  pore  pressures 
with  position  and  time  are  given  in  (a)  and  (b).  The  curve  designations 
in  (a)  refer  to  the  distance  from  the  point  of  application  of  tension 
while  in  (b)  the  designations  are  the  distance  above  the  reservoir 
water  level.  The  differences  in  the  curves  during  the  early  stages 
of  the  interval  considered,  can  be  explained  by  the  fact  that  in  the 
applied  tension  case  the  specimen  was  initially  in  equilibrium,  while 
in  the  frost  heave  tests  freezing  was  started  before  gravity  drainage 
was  complete.  With  this  exception  the  variation  of  negative  stress 
with  position  and  time  appears  to  be  similar  in  the  two  cases.  The 
similarity  of  shape  carries  through  to  the  ’’the  quantity  of  moisture 
removed  vs.  time"  and  to  the  variation  in  permeability  with  negative 
pressure,  as  shown  in  (c)-(d)  and  (e)-(f),  respectively.  In  the  frost 
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heave  tests  no  moisture  content  data  vas  obtained.  At  the 
termination  of  these  tests  some  of  the  instrumentation  was 
imbedded  in  frozen  soil  precluding  immediate  sectioning  of 
the  specimen  for  moisture  content  determinations.  Beskov  (6) 
has  shown  that  the  moisture  content  in  a  freezing  soil  specimen 
decreases  between  the  water  reservoir  and  the  frost  line.  If  the 
moisture  content  decreases  with  increasing  negative  stresses,  it 
follows  that  a  stress  distribution  such  as  that  obtained  in  the 
frost  heave  tests  should  result  in  the  rate  of  change  of  moisture 
content  increasing  in  the  direction  of  flow.  Thus,  it  is  probable 
that  there  is  similarity  in  the  moisture  content-pressure 
relationships.  The  qualitative  agreement  among  the  data  indicate 
that  the  hydrodynamics  of  flow  under  negative  pressures  is  similar 
whether  the  pressures  are  applied  or  induced. 

It  has  been  shown  by  Schofield  (29)  and  others,  that  the 
water  content  of  a  soil  is  a  continuous  function  of  the  negative 
stress  in  the  pore  water,  with  the  water  content  decreasing  as  the 
negative  stresses  increase.  The  moisture  content-tension  curves 
obtained  by  Penner  (15)  are  of  the  same  type  as  that  shown  in  Plate 
31(g).  The  moisture  content  is  most  sensitive  to  negative  pressures 
at  the  low  pressure  values.  This  indicates  that  the  moisture  content 
reduction  does  not  necessarily  involve  the  entry  of  air  into  the 
specimen.  The  same  situation  applies  to  the  pres sure -permeability 
relationships  obtained  in  the  frost  heave  tests.  This  suggests  that 
the  hydrodynamics  of  flow  under  negative  pressures  is  dependent  on 
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the  moisture -tens ion  relationships.  According  to  Gardner  (30) 
the  permeability  is  a  function  of  the  moisture  content  or  the 
moisture  tension.  This  is  apparently  accepted  by  soil  physicists 
dealing  with  problems  including  water  evaporation  from  soils. 

Due  to  the  similarities  noted  between  frost  heave  and  evaporation 
test  data,  theory  governing  the  flow  of  water  to  an  evaporation 
surface  will  be  discussed. 

"Unsaturated"  Flow  Theory 

Y/here  soil  moisture -tens ion  relationships  are  involved,  the 
term  "saturated"  is  used  to  describe  the  moisture  content  when  the 
soil  and  water  are  in  equilibrium  at  atmospheric  pressure.  Under 
pressures  less  than  atmospheric  the  moisture  content  will  be  decreased 
leaving  the  sample  "unsaturated" .  According  to  Gardner  (30)  "The  basis 
of  the  theory  of  unsaturated  flow  is  the  assumption  that  the  volume 
flux  of  water  per  unit  area  perpendicular  to  the  direction  of  flow  is 
directly  proportional  to  the  potential  gradient."  This  may  be 
expresses  as 

v-kyi  . (V.19) 

where  is  the  unsaturated  coefficient  of  permeability  which  is 
dependent  on  moisture  content  or  tension.  The  terms  v  and  i  have  the 
same  meaning  as  in  Darcy's  law  for  saturated  flow.  For  one  dimensional 
flow  in  the  z  direction 

v  =  dH/dz  .  (V.19A) 

where  dH  is  the  total  head  loss  over  a  distance  dz.  The  continuity 
may  be  derived  as  follows.  Consider  one  dimensional  flow  in  the 


z  direction  through  an  element  having  dimension  dz,  dy  and  unity 
normal  to  y-z.  If  the  superficial  velocity  of  water  entering 
and  leaving  the  element  are  v  and  v  +(3v/3z)dz,  respectively, 
then  the  rate  of  accumulation  or  depletion  of  water  in  the 
element  should  he  the  difference  multiplied  by  the  area  normal 
to  the  direction  of  flow.  i.e.  (3v/3z)dz(dy  l).  The  volume  of 
the  element  is  1  dy  dz.  Let  the  volume  of  water  in  the  element 
be  Vw  and  the  moisture  content  on  a  volume  basis  0 (i.e.  volume  of 
water  divided  by  the  volume  of  the  element).  9  =  Vw/dy  dz  and  Vw 
=  0  dy  dz.  The  volume  of  water  in  the  element  is  governed  by  flow 
to  or  from  it. 

dVw/dt  =  (<3v/3fc)dy  dz.  and  d/dt(0  dy  dz)  =  Pv/£z(dy  dz) 

Hence 


dd/dt  =  dv/dz . (V.20) 

Equation  (V.20)  expresses  the  condition  of  continuity.  This  can 
be  combined  with  equation  (V.19A)  as  follows. 


dflj/dt  =  -(d/dz)^  dH/dz) 

The  total  head  loss  (dH)  is  equal  to  the  change  in  negative  pressure 
(dh)  minus  the  difference  in  elevation  (dz). 

(V.21) 

Integrating  for  the  case  of  steady-state  flow  where  dd/dt  =  0. 

ku(dh/dz  -  l)  =  C . (v.22) 


_  a  (k,,(ah  -  dz))  _  a  (■..  ah  ) 

dS/dt  "  -  te  (  — di -  )  d?  (“»  di£  b 


The  quantity  in  the  brackets  is  the  hydraulic  gradient,  therefore  the 
constant  of  integration  C  will  be  equal  to  the  superficial  velocity  v. 
Rewriting  as  a  function  of  z  and  integrating 


-A 


dh 


1  +  v/ky 


z 


(V.23) 
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Equation  (V.23)  gives  an  unsolved  functional  relationship 
for  elevation,  negative  pressure  and  unsaturated  permeability. 

Gardner  (30)  presents  an  empirical  equation  which  relates 
unsaturated  permeability  to  negative  pressures.  The  expression 
is  as  follows: 

ku  =  a/  (h*  +  b) . (V.24) 

where  a,  b  and  n  are  constants.  He  indicates  that  the  value 
of  n  appears  to  range  between  one  and  four.  Gardner  (31)  bas 
produced  solutions  for  equation  (V.23)  using  equation  (V.2U) 
and  several  values  of  n.  The  solutions  are  as  follows: 


n  =  1  z  =  i  ln(«h  +  p)  +  K 


(V.25) 


n  =  3/2  z  = 


-  £  ln(«l^YE±h)  +  1  (V.25A) 

<*  L 5?  err*  Vhi2  w  riTTi I 


where  f-*  =  p/oC 


n  =  2  z  =  tan 


n  =  3  z  = 


o 

where  ^  J  =  P/oC 


n 


=  1+ 


=  i  [5p3Wn 


(7 

f  Vh)2 

(/=■ 

(V  p 

h  +  K 

.(  (7+h)2_') 

(/2  - 

/h  +  hc) 

(h2  + 

Phf2  +  p: 

(V.25B) 


tan" 


l(3^.?J+K(V.25C) 


+  2p372 


^h  -  pnt<i  +  p 

+PT1"l(pb/^  1)4-  v 
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where  p4  =  f3/«c 

°^=  v/a  and  (3  =#<b  +  1,  where  a  and  b  are  obtained  from  equation 
(V.24)  and  v  is  the  superficial  velocity.  These  equations  relate 
the  elevation  to  the  piezometric  pressure  and  the  superficial 
velocity. 

Unsaturated  Flow  Theory  Applied  to  Frost  Heave  Test  Data 

It  was  attempted  to  determine  a  curve  of  the  type  k^  =  a/(hn  +  b) 
which  would  fit  the  experimental  frost  heave  test  results.  Equations 
were  obtained  by  assuming  an  n-value  then  taking  and  h  values 
corresponding  to  the  coordinates  of  two  points  on  the  experimental 
curve  and  solving  simultaneously  to  determine  a  and  b.  The  derived 
equations  were  plotted  along  with  the  experimental  data  so  that  the 
closeness  of  fit  could  be  observed  (See  Plates  32A  and  32B)*  Solutions 
were  obtained  for  n-value s  from  one  to  four  because  the  flow  equation, 
(V23),  has  been  solved  for  these  values  of  n.  The  indications  are 
that  the  best  fit  curve  is  obtained  for  n  =  1.  For  n  =  1  there  was 
variation  in  the  curve  depending  on  the  experimental  points  used  in  the 
solution.  Therefore,  the  "n  =  1"  curves  are  by  no  means  an  exact 
representation  of  the  test  data.  Due  to  the  fact  that  the  fit  seemed 
to  improve  as  the  value  of  n  decreased,  n  =  0.5  was  tried,  although 
no  corresponding  solution  of  the  flow  equation  is  available.  These 
curves  (Plate  33B)  fit  the  data  rather  well  and  the  difference  in 
curves  obtained  from  two  different  pairs  of  points,  was  less  than 
that  obtained  using  other  values  of  n.  This  indicates  that  for  the 
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soil  used  and  the  conditions  imposed  in  the  present  testing  program, 
the  n-value  required  to  fit  a  curve  of  the  type  =  a/(hn  +  b)  to 
the  experimental  pressure -permeability  data  is  smaller  then  any 
considered  by  Gardner  (30)*  As  data  was  obtained  for  one  material 
only,  it  was  not  determined  whether  it  was  the  sample  material  or 
the  test  conditions  which  gave  rise  to  the  small  n-value  required 
in  the  fitting  curve.  The  lack  of  a  solution  to  the  flow  equation 
(V.23)  for  the  case  n  =  0.5>  precludes  the  continuation  of  analysis 
along  this  line;  the  derivation  being  beyond  the  scope  of  this 
dissertation. 

Analysis  of  the  Heave  Data 

The  time  and  heave  data  are  listed  in  Tables  XIII(a)  to  XIII(e). 
From  Plates  20A  to  20E  intervals  were  selected  for  which  the  rate  of 
heave  appeared  to  be  approximately  constant.  The  rates  of  heave  were 
determined  from  the  "Heave  vs.  Time"  curves  on  Plates  3^A  to  3^E. 

These  values  are  tabulated  along  with  the  corresponding  distances 
from  the  reservoir  to  the  frost  line,  in  Table  XIV.  The  rate  of 
heave  is  shown  as  a  function  of  the  length  of  flow  path,  on  Plate  35* 
For  the  first  cycle  of  freezing,  the  rate  of  heave  appears  to  decrease 
linearly  as  the  length  of  flow  path  increases.  Generally,  the  effect 
of  refreezing  was  to  increase  the  rate  of  heave.  This  effect  was  most 
obvious  in  test  12-12F.  In  12A-12D,  inclusive,  freezing  took  place 
with  the  cold  room  temperature  at  -22°F.  so  that  variations  in  the  rate 
of  heave  for  these  tests  can  be  attributed  to  effects  produced  by  the 
freezing  process  itself.  In  this  respect,  the  development  of  structure 

is  probably  involved.  Tests  12E  and  12F  showed  a  further  increase 
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in  the  rate  of  heave  but  freezing  took  place  with  higher  cold  room 
temperatures  which  could  have  provided  conditions  more  favorable 
for  ice  segregation. 

In  a  number  of  cases  the  rate  of  heave  in  a  parti cular  test 
was  constant  for  a  range  of  values  of  the  length  of  flow  path.  In 
tests  12B,  12C,  13A  and  13B  the  rates  of  heave  remained  constant 
while  the  dropping  frost  line  decreased  the  length  of  the  flow  path 
by  26  to  42$.  This  indicates  that  the  relationship  among  factors 
influencing  moisture  movements  is  such  the  net  effect  may  remain 
constant  even  though  there  is  variation  in  individual  factors. 
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TAPPING  POINT  LOCATION  DATA 


All  values  listed  in  the  above  table  are  expressed  in  centimeters. 
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TEMPERATURE  FROM  STRIP  CHARTS 
TESTS  12-12E 
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TEMPERATURES  FROM  STRIP  CHARTS 
TESTS  13. -13C 


TIME 
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40.9 

44.8 

47.8 

51.7 

55.2 

57-7 

72.0 

15.2 

23.5 

28.0 

31.2 

36.4 

40.3 

44.2 

47.2 

51.4 

55.0 

57.5 

29.8 

Test  13A 
14.2  22.1 

28.8 

33-8 

38.5 

42.3 

45.6 

48.0 

50.8 

53.6 

55.6 

38.1 

-4.0 

2.1 

13.1 

21.5 

29.I 

35.6 

40.4 

43.2 

47.8 

51.9 

54.3 

47.5 

-5.7 

-0.3 

9.9 

17.6 

26.5 

31.0 

36.9 

41.4 

46.0 

51.1 

54.1 

56.O 

-6.0 

-9.9 

9.0 

16.2 

21.9 

31.7 

35.6 

40.2 

44.7 

50.0 

53.1 

95.6 

-14.4 

-12.0 

-8.6 

-2.5 

8.4 

18.0 

25.8 

32.6 

39*5 

47.O 

51.5 

100.2 

-18.4 

-17.0 

-15.O 

-11.3 

-4.2 

4.5 

14.9 

25.9 

35.9 

45.2 

50.8 

104.1 

-18.8 

-17.5 

-15.6 

-12.2 

-5.6 

2.2 

12.1 

23.4 

32.5 

43.4 

49.7 

107.6 

-19.1 

-17.1 

-16.2 

-13.0 

-6.6 

1.2 

10.4 

21.8 

31.2 

42.0 

48.8 

108.8 

-19.O 

-18.5 

-16.0 

-13.0 

-6.6 

1.0 

10.4 

21.4 

30.8 

4l.o 

48.6 

22.1 

Test  13B 

8.0  18.0 

24.5 

31.4 

34.9 

42.6 

46.0 

47.8 

50.5 

52.2 

26.6 

5.4 

14.9 

21.0 

27.5 

32.3 

4l.O 

44.5 

46.7 

49.5 

51.4 

30.6 

1.8 

11.3 

17*9 

24.5 

29.7 

39.2 

43.1 

^5.5 

48.5 

50.6 

34.0 

1.0 

10.0 

16.3 

22.8 

27.8 

38.O 

42.1 

44.9 

47.8 

50.2 

41.5 

-0.4 

8.2 

14.4 

20.4 

25.3 

35.8 

40.5 

^3-3 

46.9 

49.4 

57.2 

-9.5 

-4.2 

3.4 

11.0 

18.0 

30.5 

36.0 

39.8 

44.5 

47.6 

100.5 

-19.7 

-18.0 

-16.8 

-13.5 

-9.8 

10.0 

19*5 

27.8 

37.8 

^3.5 

104.7 

-20.1 

-18.8 

-17.5 

-14.2 

-10.6 

8.4 

17.6 

25.8 

36.6 

43.0 

110.0 

-19.9 

-18.5 

-17.2 

-14.2 

-10.7 

7.6 

16.8 

24.4 

35.8 

42.6 

20.6 

Test 

IO" 

13C 

27.I 

35.O 

38.O 

42.0 

45.9 

47.2 

49.3 

50.6 

52.5 

53.O 

25.5 

11.4 

22.5 

30.2 

35.0 

39.4 

44.0 

45.5 

47.8 

49.5 

51.6 

52.5 

33.0 

8.0 

18.1 

26.1 

30.5 

35.7 

41.0 

42.7 

45.4 

47.3 

50.0 

51.0 

47.5 

-9.2 

-4.0 

4.6 

15.3 

24.0 

35.8 

38.0 

42.0 

44.9 

48.9 

50.3 

57*5 

-17.1 

-15.0 

-12.1 

-7.5 

-1.2 

17.2 

22.5 

30.2 

35.7 

44.9 

48.5 

80.0 

-19.0 

-17.5 

-15.6 

-12.7 

-9.5 

0.8 

11.8 

20.5 

26.4 

38.7 

43.7 

85.0 

-19.5 

-18.5 

-11 

-15.8 

-13.9 

-10.7 

-1.8 

5.2 

14.8 

3£*3, 

41.8 
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TABLE  I V  (A) 


TEMPERATURES  FROM  STRIP  CHARTS 
TEST  9 


TIME 

ELEVATION  ABOVE 

RESERVOIR 

WATER  SURFACE 

-  cm. 

hr. 

Tory 

58.0 

5?*4 

^2.9 

90.2 

4?.3 

42.7 

40.2 

35-1 

27*? 

19.9 

12.2 

1.4 

7.0 

38.0 

4l.O 

47.4 

50.7 

54.5 

57.5 

60.0 

62.0 

63.4 

65.4 

66.8 

67.6 

68.5 

68.7 

76.2 

16.0 

17.8 

21.7 

23.7 

26.4 

28.9 

31.4 

33.1 

34.5 

37.2 

42.7 

47.8 

52.4 

58.0 

90.5 

14.7 

16.0 

19.6 

21.8 

24.0 

26.1 

28.5 

30.6 

32.8 

36.4 

42.9 

48.3 

53.0 

58.9 

152.8 

13.5 

15.0 

18.8 

21.0 

23.4 

25.8 

28.0 

30.0 

32.5 

37.o 

44.1 

1*9-9 

54.3 

59.4 

I65.O 

14.2 

16.0 

19.5 

21.5 

23.6 

25.7 

27.8 

29.9 

32.1 

36.5 

43.6 

**9.3 

54.0 

60.3 

354.3 

6.4 

8.5 

13.2 

16.0 

19.5 

22.3 

25.2 

28.2 

31.5 

36.7 

44.4 

50.5 

55.3 

60.3 

388.8 

5.8 

8.2 

13.0 

15.8 

19.2 

22.0 

25.0 

27.8 

31.0 

36.6 

44.2 

50.5 

55.1 

59.8 

406.0 

5.6 

8.7 

12.2 

1^2 

18.6 

21.4 

24.5 

ZLlL 

36.1 

43.8 

50.0 

il 

NOTE:  The  temperatures  listed  were  taken  directly  from  the  strip 
charts  at  the  times  indicated  and  do  not  include  a  calibration  factor. 
Units  are  degrees  Fahrenheit. 


TABLE  IV  {&) 

TEMPERATURES  FROM  STRIP  CHARTS 
TEST  10 


TIME 

ELEVATION  ABOVE  RESERVOIR  WATER 

SURFACE  - 

cm. 

hr.  6l7cT 

97-9 

99-2 

52.8 

47.7 

T2^ 

37.6 

29.9 

22.2 

l4.6 

7.7 

30.0  15.0 

21.4 

24.2 

26.0 

31.6 

35.2 

38.3 

45.1 

50.2 

55.3 

61.3 

53.0  13.0 

18.1 

20.6 

22.0 

27.0 

31.2 

35.0 

43.8 

50.2 

56.7 

64.3 

95-5  8.8 

14.5 

17.7 

19.7 

25.5 

30.4 

34.3 

43.9 

50.5 

57.0 

64.2 

117.5  8.0 

13.4 

16.6 

18.5 

24.3 

29.0 

32.3 

41.2 

47.8 

54.0 

60.8 

145.6-10.2 

-0.5 

4.9 

8.7 

17.6 

25.2 

30.2 

38.9 

46.3 

53.0 

59.8 

153. 5  -10.4 

-1.8 

3.4 

7.0 

15.7 

23.4 

28.6 

37.7 

45.4 

52.5 

58.4 

158.0-11.2 

-2.2 

2.8 

6.1 

15.0 

22.3 

28.0 

37.0 

45.0 

51.9 

59.0 

173.4-11.4 

-3.0 

1.6 

5.0 

13.5 

20.9 

26.3 

36.2 

44.5 

52.0 

59.1 

TABLE 

iv(c) 

TEMPERATURES  FROM  STRIP  CHARTS 

TEST  11 

TIME 

ELEVATION  ABOVE  RESERVOIR  WATER 

SURFACE  - 

cm. 

hr.  61.0 

57-9 

55.3 

52.8 

47.7 

^42^ 

37.6 

29.9 

22.2 

14.6 

7.7 

29.8  19.9 

26.5 

29.2 

30.8 

35.8 

38.8 

41.8 

46.7 

51.2 

56.4 

61.4 

262.7  4.8 

14.8 

20.2 

24.3 

34.0 

39.3 

43.0 

49.0 

53.0 

57.0 

60.8 

301.5  3.7 

13.0 

18.2 

22.5 

32.0 

36.8 

40.9 

46.8 

51.0 

55.1 

58.5 

364.8  -4.5 

6.0 

12.3 

18.0 

29.8 

35.0 

39.6 

46.3 

51.1 

55.6 

59.4 

437.0-22.2 

-8.0 

0.0 

8.0 

24.2 

31.5 

37.2 

45.2 

50.7 

55.0 

58.3 

1078.8-22.3 

-10.8 

-4.8 

3.0 

16.8 

25.0 

32.2 

40.5 

46.4 

56.1 
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ELEVATION  ABOVE  FREE  WATER  LEVEL  -  Centimeters 


lit) 


ELEVATION  ABOVE  FREE  T7ATER  LEVEL  -  Centimeters 
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TEMPERATURE  PROFILES 


ELEVATION  ABOVE  FREE  WATER  LEVEL  -  Centimeters 
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ELEVATION  ABOVE  FREE  VJAiER  LEVEL  -  Centimeters 


119 


1 — 


TEMPERATURE  IROFILES 


ELEVATION  ABOVE  FLEE  WATER  LEVEL  -  Centimeters 
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VISCOSITY  -  millipoise 


121 


122 


TABLE  V 

VISCOSITY  CORRECTIONS 


TEST 

TIME 

Hr. 

ELEV. 

Cm. 

TEMP. 

°p  oc 

VISCOSITY 

Millipoise 

TOTAL  HEAD 
Cm.  of  H2O 

ELEV.  OF 

PEOST  LINE- Cm. 

(l) 

(2) 

(3) 

Gu 

(5; 

(6) 

(7) 

(0) 

(9) 

9 

76.2 

0 

60.0 

15.5 

11.3 

0 

0 

42.9 

If 

II 

12.7 

50.5 

10.3 

13.0 

1 

1 

tt 

II 

II 

17.7 

47.5 

8.6 

13.6 

1 

1 

IV 

II 

If 

22.6 

44.5 

7.0 

14.3 

21 

15 

fl 

If 

If 

26.6 

41.5 

5.3 

15.0 

27 

19 

II 

10 

173-4 

0 

63.5 

17.5 

10.8 

0 

0 

32.3 

II 

If 

13.3 

52.5 

11.2 

12.6 

24 

21 

II 

It 

If 

17.9 

47.5 

8.6 

13.6 

51 

42 

If 

fl 

II 

22.5 

42.5 

5.8 

14.8 

103 

78 

II 

It 

II 

27.3 

38.0 

3.3 

16.0 

178 

127 

If 

ft 

II 

33.1 

31.0 

-0.6 

17.9 

310 

206 

It 

11 

437.0 

0 

60.5 

15.8 

11.2 

0 

0 

40.9 

If 

It 

13.3 

54.5 

12.5 

12.2 

16 

14 

11 

It 

It 

17.9 

52.0 

11.1 

13.7 

11 

10 

11 

If 

tt 

22.5 

49.0 

9.4 

13.3 

40 

33 

II 

If 

tf 

27.3 

46.5 

8.1 

13.8 

93 

72 

II 

II 

If 

33.1 

4l.O 

5.0 

15.2 

138 

103 

II 

12B 

23.5 

0 

47.5 

8.6 

13.6 

0 

0 

12.9 

It 

ft 

2.6 

47.0 

8.3 

13.7 

9 

7 

11 

II 

ft 

5.1 

45.0 

7.2 

14.2 

35 

26 

tl 

It 

ft 

7.5 

42.5 

5.8 

14.8 

119 

109 

It 

If 

ft 

10.0 

38,5 

3.6 

15.9 

145 

126 

If 

If 

If 

12.5 

32.5 

0.3 

17.7 

379 

266 

11 

13B 

57.2 

0 

49.5 

9.7 

13.2 

0 

0 

20.8 

It 

tl 

2.6 

48.5 

9.2 

13.4 

7 

5 

If 

If 

If 

5.1 

47.5 

8.6 

13.6 

25 

18 

II 

It 

It 

7.5 

46.0 

7.8 

13.9 

86 

63 

It 

II 

II 

10.0 

44.5 

7.0 

14.3 

198 

144 

tf 

It 

II 

12.5 

42.5 

5.8 

14.8 

JZ9 

269 

If 

NOTE:  The  total  head  values  in  column  (7)  were  taken 
directly  from  Tables III(A)  toIIl(E).  The  differences  between  values 
vertically  adjacent  in  this  column  were  corrected  in  accordance  with 
equation  (II.16),  using  the  viscosity  of  water  at  20°C  as  standard. 
Starting  with  a  total  head  value  of  zero  at  the  reservoir,  the 
corrected  differences  were  added  accumulatively,  to  get  the  viscosity 
corrected  total  head  values  listed  in  column  (8).  The  average  of 
consecutive  values  in  column  (6)  was  used  as  the  average  viscosity 
for  the  corresponding  interval  along  the  specimen. 
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DETERMINATION  OF  SURCHARGE 


Weight  of  top  plate . 148  gm. 

Thrust  of  Ames  dial  shaft .  250  gm. 

Total . 398  gm. 

Diameter  of  sample  . 3.12"  =  7.94  cm. 

Cross-section  area  .  49.4  cm.2 

Fixed  surcharge .  8.1  gm./cm.^ 

Depth  of  lucite  rings . 1.00"  =  2.54  cm. 

Weight  of  lucite  rings  . 37. 5  gm. 

Weight  of  lucite  rings  per  cm.  of  length . 14.8  gm. 

Weight  of  soil  per  cm.  length  of  sample1 . 93*9  gm* 

Combined  weight  per  cm.  of  length .  109  gm. 

Combined  unit  weight  .  2.21  gm. /cm. 2 

Total  surcharge  (8  =  2.21  x  length  in  cm.)  gm./cm.^ 


Based  on  a  unit  weight  determined  using  assumed 
values  for  the  void  ratio  and  the  degree  of  saturation.  The 
maximum  error  due  to  these  assumptions  is  +  5$« 


TABLE  VII 


SURCHARGE  AT  TAPPING  POINTS 


TEST 

MANOMETER 

1 

2 

3 

4 

5 

9 

95 

116 

86 

106" 

10  73 

107 

1 16 

85 

96 

H  73 

107 

116 

85 

96 

12-12F  48 

42 

37 

31 

26 

13-13C  87 

81 

76 

71 

65 

% 


* 


r: 


•-j 
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TABLE  VIII 

VOID  RATIO  CORRECTIONS 


TEST 

TIME 

Hr. 

ELEV. 

Cm. 

EFFECTIVE 

PRESSURE 

Kg. /Cm. 2 

VOID 

RATIO 

e 

- S3 — 

1  +  e 

TOTAL  HEAD 

Cm.  of  H2O 

(l) 

(2) 

(3) 

(4) 

(5) 

(6) 

^7 r 

- IFT 

9 

76.2 

0 

0.15 

0.595 

0.132 

0 

0 

It 

II 

12.7 

0.13 

0.598 

0.134 

1 

1 

ft 

II 

17.7 

0.13 

0.598 

0.134 

1 

1 

ft 

VI 

22.6 

0.14 

0.596 

0.133 

15 

15 

ft 

It 

26.6 

0.14 

0.596 

0.133 

10 

19 

10 

173-4 

0 

0.15 

0.595 

0.132 

0 

0 

Vf 

If 

13.3 

0.15 

0.595 

0.132 

21 

21 

If 

It 

17.9 

0.18 

0.590 

0.130 

42 

42 

II 

II 

22.5 

0.22 

0.584 

0.126 

78 

77 

It 

It 

27.3 

0.29 

0.576 

0.122 

127 

123 

II 

II 

33.1 

0.42 

0.565 

0.115 

206 

194 

11 

437.0 

0 

0.15 

0.595 

0.132 

0 

0 

II 

II 

13.3 

0.l4 

0.596 

0.133 

14 

14 

II 

It 

17.9 

o.i4 

0.596 

0.133 

10 

10 

II 

tl 

22.5 

0.16 

0.593 

0.131 

33 

33 

II 

II 

27.3 

0.21 

0.585 

0.126 

72 

70 

II 

tl 

33.1 

0.24 

0.581 

0.124 

103 

99 

12B 

23.5 

0 

0.06 

O.611 

0.l4l 

0 

0 

II 

ft 

2.6 

0.06 

O.611 

o.i4i 

7 

7 

II 

tl 

5.1 

0.08 

0.608 

0.140 

26 

26 

1? 

II 

7*5 

0.16 

0.593 

0.131 

109 

105 

II 

ft 

10.0 

0.19 

0.588 

0.128 

126 

120 

It 

II 

12.5 

0.42 

0.565 

0.115 

266 

241 

13B 

57.2 

0 

0.10 

o.6o4 

0.137 

0 

0 

tl 

If 

2.6 

0.10 

0.6o4 

0.137 

5 

5 

II 

tl 

5.1 

0.11 

0.602 

0.136 

18 

18 

II 

II 

7-5 

0.17 

0.592 

0.131 

63 

62 

It 

II 

10.0 

0.28 

0.577 

0.122 

144 

137 

II 

II 

12.5 

0.46 

0.562 

0.114 

269 

244 

NOTE:  In  each  test,  the  void  ratio  at  the  reservoir  was 
taken  as  standard  for  correction  purposes.  Column  (7)  consists 
of  viscosity  corrected  total  head  values  from  Table  IV.  Differ¬ 
ences  in  consecutive  values  were  corrected  for  void  ratio  on  the 
basis  of  equation  (V.l8).  Starting  with  zero  at  the  reservoir, 
the  corrected  differences  were  added,  accumulatively,  to  get  the 
total  head  values  listed  in  Column  (8).  These  are  corrected  for 
both  viscosity  and  void  ratio  differences. 


ELEVATION  ABOVE  RESERVOIR  WATER  SURFACE 
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I 


kgxoer imental  Values 


k  Experimental  Values 


0  Values  Corrected 


©  Values  Corrected 
to  20^0 


a  Corrected  for  Viscosity 


Void  Ratio 


:perlnental  Values 

1  liA.q  fiorroi 


x  lbcperjjnental  Values 


O  Values  Corrected  to  20°C 


Corrected  to  20°C 


O  values 


&  Corrected  for  ViscoBity 


for  Vis cos lty 


Bitio 


TOTAL  HEAD 


PLATE 


wivm 


— 

L j 

WITH  AND  WITHOUT  CORRECTIONS 


% 
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HYDRAULIC  GRADIENT  AND  PERMEABILITY  COMRJTATIONS 


ELEV.  TOTAL  HEAD  HYDRAULIC  PERMEABILITY  PIEZOMETRIC 

GRADIENT  FRESSURE 

Cm.  of  H20  Cm./Sec.  x  10^  Cm.  of  H2O 


TTT 

~nr~ 

— C4T" 

~vrr~ 

~wr~ 

~rn — 

"W 

— c^j — 

Test 

11  -  437.0  hr.. 

Rate  of 

heave 

=  0.0040  Cm. /hr. 

0 

0 

0 

0.7 

0.7 

140 

140 

0 

0 

13.3 

9 

9 

1-5 

1.5 

67 

67 

22 

22 

17.9 

16 

16 

5.2 

3.7 

19 

27 

34 

34 

22.5 

40 

33 

8.3 

6.0 

12 

17 

62 

56 

27.3 

80 

62 

10 

6.4 

10 

16 

107 

89 

33.1 

138 

99 

171 

132 

Test 

12A  -  23. 

.5  hr.. 

Rate  of 

heave 

=  0.044 

Cm. /hr. 

0 

0 

0 

3.5 

2.7 

310 

410 

0 

0 

2.6 

9 

7 

10 

7.6 

110 

150 

12 

10 

5-1 

35 

26 

28 

20 

39 

55 

40 

31 

7-5 

102 

75 

49 

27 

22 

4l 

no 

83 

10.0 

225 

143 

62 

39 

18 

28 

235 

153 

12.5 

379 

241 

392 

254 

Test 

13B  -  57 

.2  hr. , 

Rate  of 

heave 

=  0.022 

Cm. /hr. 

0 

0 

0 

2.7 

1.9 

210 

290 

0 

0 

2.6 

7 

5 

7.2 

5.2 

78 

110 

10 

8 

5.1 

25 

18 

26 

18 

22 

31 

30 

23 

7.5 

86 

62 

45 

30 

12 

19 

94 

70 

10.0 

198 

137 

72 

43 

7.8 

13 

208 

147 

12.5 

379 

244 

292 

257 

NOTE: 

In  the 

pairs  of  columns  (2)- 

(3),  (4)- 

(5), 

(6)-( 7) 

and  («)-(9),  the  values  in  the  second  column  of  each  pair  contain 
void  ratio  and  viscosity  corrections,  while  those  in  the  first 
column  do  not. 
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TABLE  X 


COMPUTATION  OF  PERMEABILITY  AT  VARIOUS  TIMES  IN  TEST  13B 


ELEV. 

PIEZOMETRIC 

TOTAL 

HEAD  i 

PERMEABILITY 

PRESSURE 

HEAD 

LOSS 

Cm. 

Cm.  of  H2O 

Cm.  of  H2O 

Cm./ Sec. 

U) 

(2) 

(3J 

00  UV 

00 

0 

0 

0 

2  0.8 

2.1  x  10"5 

Time  -  22.1  Hr.  2.6 

5 

2 

3  1 

1.7  x  10-5 

Rate  of  Heave  5.1 

10 

5 

0.06l  Cm. /Hr. 

7  3 

5.6  x  10-6 

7-5 

20 

12 

26  10 

1.7  x  10"° 

10.0 

48 

38 

^3  17 

1.0  x  10'° 

12.5 

94 

81 

0 

0 

0 

5  2 

6.4  x  lO-6 

Time  -  26.6  Hr.  2.6 

8 

5 

• 

8  3 

4.2  x  10'° 

Rate  of  Heave  5*1 

18 

13 

0.046  Cm. /Hr. 

21  9 

1.4  x  10*° 

7-5 

42 

34 

54  22 

5.8  x  10-7 

10.0 

98 

88 

69  28 

4.4  x  10-7 

12.5 

170 

.157. 

0 

0 

0 

8  3 

3.1  x  lO-o 

Time  -  34.0  Hr.  2.6 

11 

8 

15  6 

1.5  x  10-fe 

Rate  of  Heave  5*1 

28 

23 

0.033  Cm. /Hr. 

46  19 

4.7  x  10-7 

7-5 

77 

69 

90  36 

2.6  x  10-7 

10.0 

169 

159 

118  47 

1.9  x  10-7 

12.5 

290 

277 
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TABLE  X  (Concluded) 


ELEV. 

PIEZOMETRIC 

TOTAL 

HEAD 

i 

PERMEABILITY 

PRESSURE 

HEAD 

LOSS 

Cm. 

Cm .  of  HpO 

Cm.  of  HpO 

Cm./ Sec. 

~tt 

T2T 

TV 

(5) 

nr — 

0 

0 

0 

11 

4 

2.4  x  10"6 

Time  -  110.0  Hr.  2.6 

14 

11 

31 

12 

8.1  x  10-7 

Rate  of  Heave  5*1 

47 

42 

0.035  Cm. /Hr. 

136 

57 

1.7  x  10-7 

7-5 

186 

178 

192 

77 

1.2  x  10-7 

10.0 

380 

370 

258 

103 

9.4  x  10-ti 

12.5 

641 

628 

NOTE:  The  rate  of  heave  was  computed  from  the  heave  curve 
on  Plate  20E.  The  total  head  values  listed  have  not  been  corrected 
for  temperature  or  void  ratio  variations.  All  the  piezometric 
pressure  and  total  head  values  are  negative. 
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TABUS  XI 

HYDRAULIC  GRADIENT  VALUES  AT  DIFFERENT  TIMES  AND  ELEVATIONS 

TEST  10. 


ELEV. 

Cm. 

TOTAL  i 

HEAD 

Cm. 

ELEV. 

Cm. 

TOTAL 

HEAD 

Cm. 

i 

(1) 

nr 

nr~ 

~Tz) — 

~JTT 

Time 

-  30*0  hr. 

Time  - 

145.6  hr. 

0 

0 

+0.3 

0 

0 

0.0 

13-3 

+4 

+1 

13.3 

0 

-1 

17.9 

+9 

+1 

17.9 

-5 

-2.4 

22.5 

+12 

0 

22.5 

-16 

-7.5 

27.3 

+13 

+1 

27.3 

-42 

-14 

33.1 

+19 

33.1 

-124 

Time  - 

95.5  hr. 

Time  - 

153.5  hr. 

0 

0 

+0.4 

0 

0 

-0.2 

13.3 

+5 

+1 

13.3 

-3 

-1 

17.9 

+8 

-1 

17.9 

-8 

-4.8 

22.5 

+4 

-2 

22.5 

-29 

-9.8 

27.3 

-5 

-4.3 

27.3 

-76 

-20 

33.1 

-30 

33.1 

-194 

Time  - 

117.5  hr. 

Time  - 

158.0  hr. 

0 

0 

40.4 

0 

0 

-0.5 

13.3 

+5 

0.0 

13.3 

-6 

-2 

17*9 

+4 

-1 

17.9 

-15 

-6.4 

22.5 

-1 

-4.4 

22.5 

-44 

-13 

27.3 

-22 

-7.2 

27.3 

-108 

-24 

33-1 

-64 

33.1 

-249 

NOTE:  Hydraulic  gradients  tending  to  produce 

upward  flow  of  water  are  shown  as  negative. 
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HYDRAULIC  GRADIENT  AFTER  TERMINATION  OF  IRRIGATION 


CAPILLARY 

SOIL  MOISTURE  CONTEIIT  CONDUCTIVITY  YIEID  TENSION 

Cm.^/Gm.  x  10^  Cm. /Hr.  Surface  Cm.  Cm.  of  Water 


<*200 


Hours 


Hour  b 


Hours 


HOTEL: 


1.  The  x  values  given  5,n  (a) 
refer  to  jtha  distance  e£  tee 
point  considered  in  the  uanpl 
from  the  point  of  applicej  t?.on 


-f  WSION  -  -Cm.  of  Water 


(After  Richard3 


2.  Trtieldfl  is  equivalent  to 


3.  Capillary  conductivity  ie 
equivalent .  to_ peroeaa  llltyv . . 


PLATE  31 


HVJinr,  ;?Rfl?  P^TTT.TS  ffTR  P'H'RLTSHHD 


A  r.(foFPART30N  OF 


RESULTS  OF  TESTS  USING  APPLIED  MOISTURE  TENSIONS 


•  T,  1 

I— — 

T  *  1 

FROSIi 

TABLE  XII 


DETERMINATION  OF  CONSTANTS  FOR  FITTING  CURVES  TO  EXPERIMENTAL  DATA 


Points  selected  from  the  experimental  curve  on  Plate  32  (a) 


1.  ^ 

= 

3.1  x  10“  6  cm. /sec 

.  h  =  6 

cm. 

2. 

= 

5.0  x  10-7  cm. /sec 

h  =  60 

cm. 

3.  K 

= 

1.7  x  10~7  cm. /sec 

.  h  =  312 

cm. 

Points 

selected 

from  the  experimental  curve  on  Plate  33  (a) 

1.  ku 

= 

2.9  x  10"  6  cm. /sec 

.  h  =  5 

cm. 

2.  k* 

= 

3.7  x  10-7  cm. /sec 

t— 

-4- 

11 

A 

• 

cm. 

3.  ky 

= 

9.4  x  10"®  cm. /sec 

• 

P* 

II 

U> 

8 

cm. 

TEST 

TIME 

POINTS  a 

b 

ku 

Hr  o 

USED 

Cm./ Sec . 

(1) 

(2) 

"T5}  (5T 

(5) 

(6) 

12  B 

23.5 

1 

-  3  56  x  10-6 

12 

56  x  10" 6/ (h  + 

12) 

12B 

23.5 

1 

-  2  32  x  lO-6 

4.2 

32  x  10-6/ (h  4 

4.2) 

13B 

57-2 

4 

-  6  28  x  10"6 

4.6 

28  x  10"  k/  (h  4 

4.6) 

13B 

57.2 

4 

-  5  18  x  10"  6 

1.2 

18  x  10"  6/ (h  + 

1.2) 

12B 

23.5 

1 

-  3  18  x  10"3 

5700 

18  x  lO“3/(h2  4 

5700) 

12B 

23.5 

1 

-  2  21  x  10" 1: 

650 

21  x  10"V(k2  + 

650) 

12B 

23.5 

1 

-  2  0.13 

4l  x  103 

0.13/ (h3  4  4l  x 

103) 

12B 

23.5 

1 

-  2  7.7 

2.5  x  106 

7*7/(h4  +  2.5  x 

106) 

12B 

23.5 

1 

-  3  27  x  10-7 

-1.6 

27  x  10-7/ (h°* 5 

-  1.6) 

12B 

23.5 

1 

-  2  32  X  10“7 

-1.4 

32  x  10-7/ (h°* 5 

-  !•») 

PERMEABILITY  -  Cm./Sec 


13* 


to 

o 

rH 

X 


^iteOMETRIC  PRESSURE 


of  Water 


PLATS 


TnTT^TWf!  fiTmvris  TQ 


\ 


TABLE  XIII  (A) 


139 


HEAVE  VALUES  FOR  TEST  9. 


TIME 

Hr. 

HEAVE 

In. 

TIME 

Hr. 

HEAVE 

In. 

TIME 

Hr. 

HEAVE 

In. 

5.0 

-0.004 

139*3 

0-974 

344.2 

1.053 

7.0 

-0.002 

140.3 

O.978 

3^5.2 

1.103 

8.5 

0.000 

141.3 

O.98I 

35^-3 

1.110 

18.3 

0.049 

143.0 

O.987 

357.3 

1.128 

19.3 

0.058 

145.0 

0.991 

360.9 

1.148 

20.3 

0.064 

146.5 

0.995 

364.1 

1.161 

21.6 

0.075 

147.6 

O.998 

366.0 

1.168 

23.5 

0.104 

149.1 

1.005 

378.8 

1.224 

24.6 

0.104 

150.0 

1.009 

381.2 

1.235 

25.6 

0.144 

151.0 

1.013 

382.6 

1.241 

26.6 

0.162 

152.0 

1.018 

384.5 

1.251 

31.7 

0.245 

152.8 

1.022 

386.0 

1.258 

44.8 

O.436 

162.3 

1.081 

387.2 

I.265 

47.2 

0.464 

163.8 

1.090 

388.8 

1.272 

69.5 

0.690 

165.0 

I.O96 

390.8 

1.280 

76.2 

0.755 

320.0 

0.914 

392.4 

1.295 

90.5 

0.859 

330.3 

0.979 

403.0 

1.348 

92.3 

0.872 

331*1 

0.984 

406.0 

1.363 

93.6 

0.881 

332.0 

O.989 

408.3 

1.375 

96.2 

0.899 

333*0 

0.995 

410.3 

1.383 

98.2 

0.912 

335*1 

1.008 

412.0 

1.390 

102.8 

O.929 

336.2 

1.013 

413.1 

1.395 

114.5 

O.938 

337.1 

1.018 

4l6.o 

1.407 

117.1 

0.934 

338.2 

1.023 

426.3 

1.465 

120.0 

0.932 

339.0 

1.029 

428.7 

1.478 

123.4 

O.932 

340.1 

1.036 

432.3 

1.496 

125.0 

0.933 

341.5 

1.042 

434.3 

1.503 

126.2 

0.935 

342.5 

1.046 

438.0 

1.515 

127.2 

0.939 

450.3 

1.574 

128.0 

0.942 

452.2 

1.584 

138.3 

O.97O 

454.3 

1.595 
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table  XIII  (B) 


HEAVE  VALUES  FOR  TEST  10. 


TIME 

Hr. 

HEAVE 

In. 

TIME 

Hr. 

HEAVE 

In. 

6.0 

0.080 

112.6 

0.234 

16.5 

0.150 

113.7 

0.234 

19.0 

O.I63 

115.6 

0.235 

21.6 

0.175 

117.5 

0.235 

24.0 

0.179 

119.4 

0.235 

26.4 

0.180 

121.4 

0.235 

30.0 

0.191 

126.2 

0.235 

40.4 

0.216 

136.4 

0.235 

42.5 

0.217 

136.4 

O.236 

46.0 

0.217 

141.3 

0.236 

50.5 

0.218 

143.4 

0.236 

54.7 

0.219 

145.6 

0.238 

67.3 

0.221 

148.9 

0.249 

71.0 

0.221 

153.5 

O.278 

74.9 

0.222 

160.5 

0.342 

78.9 

0.223 

163.0 

O.367 

88.5 

0.224 

165.5 

0.390 

90.5 

0.224 

167.9 

o.4n 

93.5 

0.224 

169.2 

0.421 

95.5 

0.225 

171.7 

0.443 

97.4 

0.225 

173-4 

0.456 

99.6 

0.226 

17^.3 

0.463 

101.7 

0.228 

.  . 
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TABLE  XIII  (C) 
HEAVE  VALUES  FOR  TEST  11. 


TIME 

HEAVE 

TIME 

HEAVE 

TIME 

HEAVE 

Hr. 

In. 

Hr. 

In. 

.... 

In. 

7.0 

0.130 

275.0 

0.823 

718.7 

1.093 

11.0 

0.211 

287.5 

0.852 

729.2 

1.120 

14.6 

0.249 

292.0 

O.858 

742.8 

1.140 

22.7 

0.341 

297.6 

0.866 

752.6 

1.151 

27.7 

O.368 

301.3 

0.869 

769.4 

1.180 

29.8 

O.378 

310.5 

0.880 

779-4 

1.201 

33.8 

0.401 

316.1 

0.884 

790.5 

1.216 

46.8 

0.454 

334.7 

0.902 

891.2 

1.256 

49.7 

0.455 

340.8 

0.913 

900.7 

1.257 

53.3 

0.453 

345.5 

0.921 

910.6 

1.266 

56.7 

0.472 

349.5 

O.929 

916.3 

1.272 

60.6 

0.490 

358.9 

0.938 

921.3 

1.274 

71.O 

0.553 

364.8 

0.942 

934.7 

1.308 

75.8 

O.567 

369.6 

0.946 

939.9 

1.318 

80.6 

0.575 

382.7 

0.959 

945.8 

1.325 

84.0 

O.58I 

390.7 

O.969 

958.4 

1.343 

94.7 

0.619 

397-0 

0.975 

963.1 

1.349 

99.8 

0.625 

406.6 

0.989 

972.2 

1.359 

102.7 

0.626 

417-3 

0.997 

982.6 

1.372 

17O.O 

0.480 

430.5 

1.013 

987.2 

1.376 

175.3 

0.473 

437.0 

1.023 

996.6 

1.386 

181.1 

0.484 

445.8 

1.036 

1006.6 

1.395 

191.1 

0.557 

455.0 

1.033 

1011.1 

1.399 

195.9 

0.579 

462.3 

1.034 

1014.9 

1.402 

201.6 

0.599 

583.0 

O.816 

1020.4 

1.406 

214.8 

0.675 

598.5 

0.888 

1030.6 

1.417 

219.7 

0.696 

612.5 

0.946 

1035.3 

1.423 

223.7 

0.706 

623.7 

0.983 

1043.6 

1.431 

227.1 

0.714 

633.2 

1.012 

1054.6 

1.441 

238.5 

0.747 

649-7 

1.046 

1066.0 

1.452 

244.8 

0.759 

670.5 

1.056 

IO78.7 

1.468 

250.4 

0.761 

680.5 

1.059 

1091.3 

1.480 

262.7 

0.797 

094.7 

1.085 

1127.1 

1.515 

265.8 

0.806 

700.0 

1.093 

1155-4 

1.538 

270.7 

O.813 

r  jL'f 
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TABLE  XIII  (D) 
HEAVE  VALUES  FOR  TESTS  12-12F 


TIME 

HEAVE 

TIME 

HEAVE 

TIME 

HEAVE 

Hr. 

In. 

In. 

In. 

Test  12 

Test  12A 

Test  12D 

2.2 

-o.oo4 

5.5 

0.000 

”5-8 

0.006 

5-0 

-0.005 

7-5 

0.051 

6.5 

0.026 

17-4 

0.000 

8.5 

0.085 

7*3 

0.051 

19.4 

0.001 

10.2 

0.145 

9-1 

0.095 

22.9 

0.001 

11.3 

0.198 

10.6 

0.195 

25.3 

0.001 

11.8 

0.220 

11.0 

0.217 

29.3 

0.007 

14.4 

0.291 

12.0 

0.281 

41.5 

0.187 

23.4 

0.450 

13.0 

0.342 

42.8 

0.206 

13.5 

O.366 

44.5 

0.222 

Test  12B 

23.0 

O.696 

46.1 

0.218 

3.0 

0.000 

23.5 

0.706 

47.5 

0.189 

5.0 

0.011 

49.0 

0.147 

6.0 

0.033 

Test  12E 

55-0 

0.101 

7-2 

0.065 

"XT' 

0.003 

65.5 

0.097 

8.1 

O.O98 

7.8 

0.058 

67.3 

0.096 

10.4 

0.152 

9.6 

0.128 

70.0 

0.095 

12.0 

0.249 

20.6 

0.519 

72.0 

0.096 

13.5 

0.280 

21.4 

0.544 

74.0 

0.099 

23.5 

0.458 

22.5 

O.569 

77.3 

0.115 

24.0 

0.469 

23.4 

0.595 

79.1 

0.155 

89.5 

0.353 

Test  12C 

Test  12 F 

90.0 

0.359 

3.7 

0.000 

7.0 

0.007 

91.3 

0.361 

6.4 

0.021 

17.5 

0.254 

92.5 

0.357 

12.6 

0.289 

18.3 

O.276 

94.1 

0.367 

19.4 

0.407 

19.5 

0.319 

95-1 

0.4l6 

20.0 

0.4l6 

20.0 

0.332 

96.0 

0.432 

21.0 

0.434 

20.5 

0.349 

99-2 

0.483 

22.0 

0.448 

20.8 

0.356 

100.2 

0.498 

23.0 

0.464 

22.3 

0.394 

113.5 

0.686 

24.0 

0.477 

23.0 

0.412 

114.5 

0.702 

115.5 

0.719 
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'CABLE  X.II J  (E) 

HEAVE  FOR  TESTS  13-13C 


TIME 

Hr. 

HEAVE 

In. 

TIME 

Hr. 

HEAVE 

In. 

TIME 

Hr. 

HEAVE 

In. 

Test  13 

48.5 

0.703 

36.0 

O.632 

8.8 

0.001 

49.6 

0.717 

46.0 

O.709 

20.5 

0.074 

50.5 

0.728 

V7.3 

O.723 

21.7 

0.138 

52.5 

0.757 

49.4 

0.742 

23.1 

O.169 

53.5 

0.774 

52.0 

•  O.769 

24.5 

O.I98 

54.5 

O.788 

54.5 

O.787 

25.7 

0.226 

55.5 

0.800 

57.2 

0.813 

28.0 

O.27O 

59-7 

0.845 

70.0 

0.891 

28.8 

0.282 

62.1 

0.870 

76.3 

0.964 

44.5 

0.508 

71.5 

0.931 

78.5 

O.98I 

47.2 

0.536 

72.5 

O.938 

83.0 

1.021 

49.2 

0.553 

73-5 

0.944 

97-3 

1.130 

51.5 

0.57^ 

74.5 

0.950 

99.0 

1.151 

53.0 

O.585 

76.2 

O.967 

100.5 

1.173 

55-2 

0.606 

77-5 

0.975 

110.0 

1.312 

57-2 

0.622 

79-8 

0.976 

118.2 

1.426 

68.5 

O.705 

95-8 

0.977 

69.8 

0.712 

98.1 

1.009 

Test  13C 

72.0 

0.721 

100.2 

1.045 

20.6 

0.6l0 

73-2 

0.722 

101.4 

1.066 

22.0 

O.634 

74.6 

O.73O 

102.5 

1.087 

23.6 

0.668 

107.6 

1.295 

25.5 

0.706 

Test  13A 

108.0 

1.300 

27.5 

0.736 

6.4 

0.017 

108.8 

1.311 

29.1 

0.755 

8.2 

O.O72 

33.0 

0.784 

10.1 

0.115 

Test  13B 

44.6 

0.936 

12.6 

0.205 

6.5 

0.022 

47.5 

0.975 

23.5 

0.446 

10.0 

0.140 

53.0 

1.001 

25.4 

0.464 

12.0 

0.162 

56.5 

1.025 

26.2 

0.473 

22.1 

0.421 

68.5 

1.045 

27.O 

0.483 

23.1 

0.442 

73.2 

1.120 

28.2 

0.496 

24.0 

0.460 

75.1 

1.140 

29.1 

0.504 

25.2 

O.476 

76.2 

1.146 

29.8 

0.510 

26.6 

0.505 

79-0 

1.187 

32.0 

O.528 

27.7 

0.520 

81.0 

1.210 

34.3 

0.554 

29.0 

0.545 

82.0 

1.222 

36.1 

0.578 

30.6 

O.567 

83.0 

1.234 

38.1 

0.603 

34.0 

0.609 

84.0 

1.252 

47.5 

0.693 

35.0 

0.622 

85.O 

I.269 

i 


T 


HEAVE  -  Inches 


141 


HEAVE  -  Inches 


145 


14(i 


Inches 


147 


TABLE  XIV 


DISTANCE  FROM  WATER  RESERVOIR  TO  FROST  LINE 
FOR  VARIOUS  RATES  OF  HEAVE 


TEST 

TIME 

hr. 

DEPTH  OF 
FROST  LINE 
in. 

RESERVOIR 

TO  FROST  LINE 
in. 

RATE  OF 
HEAVE 
cm. /hr. 

9 

120-217 

8. 4-8.9 

16.0-15.5 

0.009 

It 

1+00-459 

9. 4-9. 6 

15-0-14.8 

0.013 

tr 

800-885 

10.2-11.4 

14.2-13.0 

0.011 

10 

80-148 

8.8-11.0 

15.7-13.5 

0.001 

ir 

148-180 

11.0-12.0 

13*5-12.5 

0.022* 

11 

190-288 

4. 2-5. 9 

20.3-18.6 

0.005 

u 

620-690 

7.8-8. 1 

16.7-16.4 

0.002 

ii 

1100-1188 

10.2-10.9 

14.3-13.6 

0.002 

12 

97-117 

1.0-1. 4 

6. 6-7.0 

0.032 

12A 

23.28 

2.9 

5-1 

O.O36 

12B 

12.25 

1. 2-3.O 

6. 8-5.O 

0.044 

12  C 

13-24 

1. 5-3*3 

6. 5-4. 7 

0.04l 

12D 

23-25 

3*1 

4.9 

0.047 

12E 

21-28 

2. 5-3*0 

5. 5-5-0 

0.070 

12  F 

22-25 

3*5 

4.5 

0.06l 

13 

45-65 

2. 6-3. 8 

12.4-11.2 

0.021 

13A 

25-75 

2. 7-6.0 

12.3-9.0 

0.026 

13B 

35-100 

5* 5-9*5 

9. 5-5. 5 

0.022 

13C 

50 

5-4 

9.6 

0.020 

*  The  frost  line  had  penetrated  the  silt. 
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CHAPTER  VI 


CONCLUSIONS  AID  RECOMMENDATIONS 

Conclusions 

1.  The  fluctuations  of  the  heave,  depth  of  frost  line  and 
piezometric  pressure  values  that  were  observed  in  the  frost  heave 
tests,  appear  to  support  Martin's  "rhythmic  ice  banding"  hypothesis. 

2.  The  magnitude  of  the  fluctuations  tended  to  increase  as 
the  rate  of  penetration  of  the  frost  line  decreased,  indicating  an 
increasing  instability  in  factors  controlling  the  moisture  flow. 

3.  An  "overall  steady-state"  condition  appeared  to  exist 
where,  averaged  over  an  extended  period  of  time,  the  fluctuating 
quantities  could  be  defined  by  increasing  linear  functions  of  time. 

4.  The  piezometric  pressures  exceeded  the  capillarity  of  the 
material  in  the  late  stages  of  one  test  and  in  the  vicinity  of  the 
frost  line.  As  these  data  were  excluded  from  the  pressure  analysis, 
the  flow  characteristics  obtained  could  not  be  attributed  to  the 
entry  of  air  into  the  specimen. 

5.  It  was  found  that  after  taking  into  account  variations  in 
the  viscosity  and  void  ratio  along  the  flow  path,  the  characteristics 
of  the  moisture  movements  in  the  frost  heave  tests  could  not  be 
explained  in  terms  of  Darcy  type  flow,  despite  the  fact  that  the  degree 
of  saturation  remained  constant  at  96$* 

6.  A  comparison  of  the  frost  heave  test  results  with  published 
data  indicated  a  similarity  in  the  hydrodynamics  whether  the  water 
was  removed  from  the  specimen  by  an  ice  lens,  evaporation,  or  by  an 
applied  negative  pore  pressure. 
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7 •  The  fact  that  the  experimental  pres sure -permeability  data 
can  be  fitted  by  a  curve  of  the  type  k^  =  a/(hn  +  b)  indicates  that 
the  unsaturated"  flow  theory  may  apply  to  the  moisture  movements  that 
occur  during  ice  segregation. 

8.  There  are  insufficient  data  to  indicate  whether  the  small 
n- value  that  produced  the  best  fitting  curve  is  due  to  the  specimen 
material  or  the  conditions  imposed  during  the  frost  heave  tests. 

9*  For  the  material  tested,  first  cycle  freezing  indicated 
that  the  rate  of  heave  was  inversely  proportional  to  the  distance 
between  the  water  reservoir  and  the  frost  line. 

10.  The  rate  of  heave  tended  in  increase  with  each  cycle  of 
freezing.  No  maximum  was  established  at  the  end  of  seven  cycles. 

Recommendations 

Further  investigation  should  be  carried  out  with  respect  to  the 
application  of  "unsaturated"  flow  theory  to  the  frost  heave  problem. 
Frost  heave  tests  on  a  variety  of  materials  would  help  to  establish 
the  range  of  n-values  that  apply  to  the  pressure-permeability  curves  for 
flow  due  to  ice  segregation.  The  possibility  of  obtaining  solutions 
to  the  flow  equation  (V.23),  other  than  those  given  by  Gardner  (31), 
should  be  investigated.  It  would,  be  possible  to  check  the  validity 
of  "unsaturated"  flow  theory  for  the  case  of  moisture  movements  due  to 
ice  segregation,  once  solutions  to  the  flow  equation  were  available 
corresponding  to  n-values  determined  from  the  experimental  pressure- 
permeability  data. 

It  is  suggested  that  the  influence  of  temperature  on  the  flow 


characteristics  can  be  determined  by  comparing  data  from  an  evaporation 


test  with  that  from  a  frost  heave  test  on  the  same  material.  It 
is  possible  that  such  tests  could  be  carried  out  concurrently  on  the 
same  specimen. 

As  the  instrumentation  required  makes  it  difficult  to  compact 
the  test  specimens,  it  is  suggested  that  consolidation  effects 
during  the  test  could  be  reduced  by  letting  the  sample  consolidate 
under  a  surcharge  prior  to  the  application  of  freezing  temperatures. 
For  non- swelling  materials,  a  surcharge  of  the  same  order  of  magnitude 
as  the  maximum  effective  stresses  anticipated  should  eliminate  most 
of  the  void  ratio  variation. 

In  frost  heave  testing  it  would  be  desirable  to  measure  the 
rate  at  which  water  enters  the  specimen.  By  comparing  this  quantity 
with  the  rate  of  heave  it  could  be  determined  if  the  flow  as  steady- 
state  or  not. 

The  apparatus  should  be  modified  to  permit  the  removal  of  the 
specimen  at  the  end  of  the  test.  This  would  allow  the  frozen  portion 
to  be  inspected  for  the  location  of  ice  lenses.  The  unfrozen  portion 
could  be  immediately  sectioned  to  determine  the  moisture  contents 
corresponding  to  the  terminating  pressure  conditions. 

It  is  suggested  that  the  possibility  of  using  some  electrical 
resistance  method  to  determine  the  moisture  content  in  the  test 
specimen  during  the  test  be  investigated.  By  obtaining  moisture 
content  measurements  at  the  manometer  tapping  points,  the  moisture- 
tension  relationships  could  be  obtained  from  the  frost  heave  tests. 

Frost  heave  tests  with  pressure  measurements  should  be  tried 
using  various  additives  to  see  if  changes  in  such  characteristics  as 
viscosity  and  surface  tension  are  reflected  in  the  pressure  readings. 
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DETAIL  DRAWING  No  2 
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APPENDIX  "B" 


UNIVERSITY  of  ALBERTA 

DEP’T  of  CIVIL  ENGINEERING 

SOIL  MECHANICS  LABORATORY 

ATTERBERG  LIMITS 

PROJECT  Ifil 

SITE 

SAMPLE  SILT  Cz)  * 

■LOCATION  : 

HOLE  DEPTH 

TECHNICIAN  DATE 

Liquid  Limil 

iol  No- 

1 

2 

3 

4 

C 

o.  of  Blows 

27 

23 

25 

14 

12 

11 

ontainer  No- 

V26 

V22 

V5 

V3 

V21 

t-  Sample  Wet  t  Tare 

82.8514 

81.6669 

71.4465 

75.7723 

81.1934 

*6.5025 

t-  Sample  Dry  +  Tare 

30.9972 

79.6965 

70.2098 

73.2513 

73.3341 

83.  5720 

t-  Water 

3 

to 

• 

r-l 

1.9704 

1.2367 

2.5210 

2.3593 

2.9305 

are  Container 

73.5156 

71.9736 

65.2289 

64.0741 

70.3142 

73.0865 

t-  of  Dry  Soil 

7.4316 

7.7229 

4.9809 

9.1772 

8.5199 

10.4855 

oist ure  Content  ii j°/o 

_  24.3 

.  25.5 

.  24.8 

.  .27,4  . 

 27.7  

23.0 

|W  1  -  Sample  Wet+Tare  33.6706 


P I  aslic  Li  mi  1 


■Irjq.L.No. 


Container  No- 


W4  Sample  Dry+Tare  38.4592 


Wt  Water 


Tore  Container  37.4632 


Wt  of  Dry  Soil 


Moisture  Content 


3S«?235| 


0.2114 


0.9910 


21.3 


33.7199 


0.9688 


21.0 


0.2036  0.1 


37.7911  35.680; 


p.ra; 


22.: 


Shrinkage  Limit 


Trial  No 


Container  No 


Wt  Sample  Wct+Tore 


V/t  Sample  Dry  -t-Tarc 
MH£j i 


Tore  Container 


Wt-  of  Dry  Soil  W< 


Moisture  Content  -ur7< 


.Vo I  •  Container 


j*  Vol-  Dry  3 o i!  Pat  Va 


Shrinkoqe  Vol-  V-Vo 


Shrinkage  Limit 


**  -  *  ,0°) 


Description  of  Sample  _ 

Slightly  pi-. stic  silt.  According  to 

the  Airfield  Classification  Sy3tec 

.  is  ML. 


Remorks  The  sample  had  been  air 
dried  prior  to  being  ~;3ed  ir.  V:  is 

est.  
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UNIVERSITY  of  ALBERTA 
DEP’T-  of  CIVIL  ENGINEERING 
SOIL  MECHANICS  LABORATORY 


bate 

Temp- 

0 

C. 

Time 

Elapsed 

Time 

r; 

17 

O 

3 

D 

m.m. 

Rh+mrcd 

w  % 

umi  t 

W  % 
Basis  Orig 
Sample 

Remarks 

lee.  11 

0920 

30s 

23.8 

29.1 

0.053 

29.1 

68.0 

1  " 

lm 

27.5 

27.3 

0.042 

27.3 

64.0 

1  " 

2m 

25.3 

25.6 

0.030 

25.6 

. 

1  " 

21.5 

4m 

23.9 

24.2 

0.021 

24.2 

56.6 

1  n 

-21,4  . 

8m 

17.0 

17.3 

0.016 

17.3 

40.5 

1  n 

M 

15m 

14.0 

14.3 

13 

14.3 

33.4 

1  11 

N 

0950 

30m 

9.9 

10.2 

0.0091 

10.2 

23.8 

1  " 

21.2 

1020 

lh 

7.2 

7.5 

0.U065 

7.4 

17.3 

1  " 

21.2 

1120 

2h 

 5.3 

5.6 

0.0048 

12.9 

1  ” 

20.5 

1320 

4 il 

4.3 

5.1 

.  0034 

4.9 

11.5 

fee.  12 

19.7 

0930 

24h 

3.2 

3.5 

0.0014 

2.8 

6.5 

be,  13 

19.5 

0830 

47h 

2.9 

3.2 

0.0010 

2.5 

5.9 

PROJECT 
SITE 


5AM  P  LE  SILT  T2T 


LOCATION 

HOLE 


DEPTH 


W5 


ydrometer  Nos 

'  <*>  =  "wf-  •  mt~  cd^’"  — 


and 

0*7 


Graduate  No. 


jeniscus  correction  =  cm  = 


0.3 


(Rh-rmt-  cd) 
and 


lisp er sing  agent  used  6f  Calyon  solu  ion 


_ respectively 

Amount  10  cc. 


orrection  for  change  in  density  of  liquid  due  to  addition  of  dispersing  agent  =  cd 


cd  = 


0.3 


pecific  Gravity  of  Solids  =  Gs 


and 

’?0 


jespe  ctively 


le^rrintinn  nf  ^rimnls 

Method  of  Preparation 

Isli.yhtly  Dlastic  silt.  According  to 

Treated  with  Hydrogen  Peroxide  and 

1  the  Airfield  Classification  System 

warmed  in  oven.  Soaked  with  distilled 

lit  is  an  ML. 

water  over  nirht.  then  added  Caldoji, 

Remarks  Dry  wfeirht  of  the  sample _ 

used  was  computed  on  the  basis  of  the 

aterial  recovered  after  the  test. 

]  Initial  Moisture  Content 

Dry  Weight  of  Sample 

1  Inn  t  n  i  n  p  r  Mn 

Container  No 

It.  Sample  Wet  +  Tare - — - 

Wt.  Sample  (W4v(/Dry)  t  Tare  — ■-•259,9  BO* — 
Tare 

lit  Water 

Wt-  /  Dry)  Soil  _ _L _ 

lire  Container 
lit-  of  Drv  ^nil 

Dry  Weight  from  Initial 

100  *  Wt- Wet  Soil  _ 

llitial  Moistur©  ¥o  . - - - 

MoiSTure-  ,00  +  Init- Moisf  9b 

PLATE  3B 

. 


niliol  Dry  Weight 
etained  No.  4 
ore  No - 


Vt-  Dry  +  Tare. 

are _ 

Iff-  Dry _ 


UNIVERSITY  of  ALBERTA 

DEP'T-  of  CIVIL  ENGINEERING 
SOIL  MECHANICS  LABORATORY 

SIEVE  ANALYSIS 


otal  Dry  Weight 
Sample 


Sieve 

No. 


Size  of  Opening 


Inches 


Jt 


* 


*8 


•185 


Mm. 


19-10 


9-52 


4-76 


PROJECT  _  .  ^ 

SITE 

ion 

sample 

CALLAR7  GRAVE L 

LOCATION 

HOLE 

DEPTH 

TECHNICIAN 

DATE 

Weight 

Retained 

gms. 


8824 


2320 


2790 


2220 


3031 

"3151" 


ft. 


Th  on 
gms. 


26  589 


24  269 


21  479 


19  259 


16  228 


13  067 


/•  cr  i  iH/ii 


100.0 


.ZSii— L- 

68.6 


60,7 
54*  4~ 


45.9 


36,9 


Passing 


13067 


nitial  Dry  Weight 
assing  No- 4 
are  No- _ 


i*iI4. 


079 


2  000 


293.6 


1  260 


100.0 

31.1 


20 


0331 


•840 


.423^ 


36 


53.7 


Vt-  Dry  +  Tare. 

are _ 

|Vt-  Dry _ 


40 


0165 


6  0 


•0097 


14.20. 

•250 


-522*2. 


307 


19.8 


173.9 


133 


8.6 


100 


•0059 


149 


52.0 


81 


5.2 


200 


•0029 


•074 


27.7 


53 


3.4 


36.9 


28.1 _ 

19.3 
7.1 


2.1 

1.8 


1.2 


Passing 


200 


53.4 


Description  of  Sample  According  to  the 
Airfield  Classification  System,  this 


material  would  be  called  GF.  The 


PUBLIC  ROADS  Classification  would 


probably  be  A  -  2. 


Time-  of  Sieving 


20  minutes 


Method  of  Preparation 


Remarks 


Gravel  Sizes  Sand  Sizes  Is* 

3"  2" 

j 

"  2 

Sieve  Size  * 

*4  *10  2 

A  Jt 

Q_ 40  60  100  200 

n — 

T~ 

I 

Hie 

: 

. 

Do  ^*3  rrvrr 

v 

V 

— 

^ - — - - 

* 

'C  0T 

rv s\~  m  t 

_ 

50 

K 

— 

— ^ 

_ 

L 

 . 

_ 

- - 

. 

— 

r.~  ,J ...  -  1 

— 

> 

. 

— 

- - 

.  — —  •  - 

r  rx  .. 

_ L_ _ l _ X- 

7 

T 

n  r 

L  ,  i  — 

JILi-LD 

t—l—i _ 

M 

iL 

7- 

M 

>  i 

itl  ; 

LJ — 1 - * — 

Mi  i  i  •  H  i 

n 

— 

±UL 

i 

iOC 


so 

30 

70 

60 

&50 


30 


20 


10 


Grain  Size- Mm. 


Note;  M  IT  Grain  Sizo  Scale 


PLATS  40 


PROJECT _ 

SITE 

ToT 

-SAMPLE  _ 

GRAVEL  (T) 

LOCATION 

HOLE 

DEPTH 

TECHNICIAN 

DA'i  F. 

UNIVERSITY  of  ALBERTA 

DEP’T-  of  CIVIL  ENGINEERING 

SOIL  MECHANICS  LABORATORY 

SIEVE  ANALYSIS 


otol  Dry  Weight 

if  Sample 


nitial  Dry  Weight 
detained  No.  4 
are  No _ 


N\-  Dry  +  Tare, 
are _ 


Vt-  Dry. 


Passing 


Sieve 

No. 


4 


Size  of  Opening 


Inches 


1- 


3/8 


185 


M  m. 


1910 


9-52 


4  76 


Retained 

gms. 


340 


201 


jssl 


403 


3  101 


‘o  t  a  I  Wt 
mei  Thon 
gm*. 


4  $14 


4  174 


3  333 


3  101 


Percent 

Finer  Thon 


100.0 


Vo  Finer  Than 
_/ig. 
Sample 


92  .T 


86.1 


77.7 


68.7 

5°»o 

i3.2_ 


68.7 


47.1 


nitial  Dry  Wei  ght 
Massing  No-4 
are  No  _ 


l  0 


•079 


2  000 


These 


values  'nor* ;  based  on 


20 


0331 


•840 


Wt-  Dry  +  Tare 

Tare _ 

Wt-  Dry _ 


the  grad i  nr  found  for  the  mlmis 


40 


•0165 


li20 


t/ 4  material  in  the  CALT-AP.Y 


6  0 


■0097 


■250 


GRAVEL 


100 


•0059 


149 


36.2 

25.3 
_E.l. 
_4?l. 

2,5 


200 


•0029 


•074 


lA 


Possing 


200 


Description  of  So m pie  According  to  the 
Airfield  Classification  System,  this 

material  would  be  called  GF.  The 


Public  Roads  Classifies  "ion  r.-ould 

probably  be  A  -  2, 


Time  of  Sieving  20  minutes 


Method  of  Preparation 


Remarks  This  material  is  the  sar.e  as 
CALGARY  GPJIVEL  G)  ,  with  an  excess 
of  material  passing  the  #4  sieve. 


Gravel  Sizes 


Send  Siz^s 


i 


(00 


90 

80 

70 

c 
o 

f  60 


eve 


Si^e 


T 

r  ' 

i  ii  1. 1 

TTTTl 

1.  - 

11 

1 

J 

h  ~  f 

L 

i  ii  1 

L 

DIE 

JZL 

.  . 

.  .  _ 

D10 

060‘- 

CU  = 

0.44  nrvm 

i 

I 

! 

-  .  -f- -■  - -  .4 

13  m  rr 

. 

I A 

_ 1 _  -4 

: 

lv 

1 

1 

_ 

l_^_— 

— 

- 

i 

L. 

_ 

_ 

_ 

... _ i 

\  i 

| 

X 

• 

: 

Py 

— 

-T\ 

L 

_ 

n  Li  n 

^LL-L . 

J-L 

r 

i 

Li  1 

i 

LL 

ii  ti  1 

LjJl3^±==: 

±h 

■LL 

4.50 
c 

[I 

40 

c 
1! 


030 
20 
10 
0 


Grain  Size-Mrrv. 


Note:  M  IT  Grain  Size  Scale 


PLATE  4-1 


t  ■■ 


h. 


UNIVERSITY  of  ALBERTA 

SITE 

Thh 

DEP'T  of  CIVIL  ENGINEERING 

SAMPLE 

Silt 

SOIL  MECHANICS  LABORATORY 

LOCATION 

CONSOLIDATION  DATA 

HOLE 

DEPTH 

TECHNICIAN 

DATE 20  ,I»o.59 

3inq  Doto 


ing  No  _ 

Weight  gms. 


C1 


■2Q?tQ 


hickness  ins  _ 
tameter  ins.. 
Area  sq  cm's 

Machine  Data 


2.60 


34*2 


Machine  No-. 


Multiplication  Factor _ 

VI  Block  +  Stone  +  Ball 


100 


gms. 


421 


Consolidation  Sample  WeiQhts 

Wt. Tore  ♦  Ring  +  Soil  -rWoter(End)  gms  1109. 1 

Wt  Tare  r Rmg +•  Soil  lEnd)  gms  _ 

Wt-  Tare  (Tare  No  V62  )gms _ 

Wt  Ring  +  Soil  +■  Woler  (End)  gms - 


Wt  Ring  +  Soil  ♦  Woter  (Stort)gms_ 

Wl  Ring  +S01I  gms _ 

Wt.  Soil  gms  _ 

Woter  (End)  -  _28*6 _ gm< . :  19.9 

Woter  (Stort)  .  3.6*7 _ gms  -.23. m3. 


10SC. 5 

1075*1. 

1083.6 


1046.9 


U3.9  . 


)escription  of  Sample 


ML  -  See  Atterherg  Limits 


)ate,Time 
_oad,Dial 
emarks 


Date,  Time 

Date,  Time 

Date ,  Time 

D  ate 

,  Time 

Date,  Time 

Load,  Dial 

Load, Dial 

Load, Dial 

Load 

,  Dial 

Load,  Dial 

Remarks 

Remarks 

Remarks 

Remarks 

Remarks 

20  Jan.  59 

lm  .7363 

21  Jan.  59 

2ra 

.7080 

21  Jan.  59 

20  gm. 

2m  46 

20  gm. 

4m 

80 

200  gm. 

15:05 

4m  28 

10:45 

9m 

30 

15l47 

18s  .7848 

8m  13 

6s  .7104 

21  Jan.  59 

63  .7062 

30s  40 

15m  06 

12s  06 

400  gm. 

12b  61 

Ira  22 

30m  00 

30s  10 

11:15 

30s  60 

2d  01 

50m  .7298 

lm  13 

6e 

.7060 

2m  60 

Added  water 

21  Jan.  59 

2ra  17 

12b 

58 

4m  59 

4m  .7760 

200  go. 

4m  19 

30s 

54 

9m  59 

8m  06 

10:10 

Sm  20 

lm 

51 

21  Jan.  59 

15m  .7631 

12s  .7248 

21  Jan.  59 

2m 

50 

400  gm. 

30m  .7544 

30s  35 

50  gm. 

4m 

49 

15:58 

75m  24 

lm  24 

10:54 

8m 

46 

6s  .7033 

Covered  wit' 

1  2m  11 

63  . 711 5 

15m 

44 

12s  36 

water. 

4m  •  00 

12s  15 

3  On 

43 

30s  34 

120m  .7519 

8m  .7192 

In  13 

142m 

22- 

lm  32 

1060m  13 

15n  88 

2m  11 

21  Jan.  59 

2m  30 

21  Jan.  59 

21  Jan.  59 

4m  11 

20 

gm. 

4m  30 

50  gm. 

400  gm. 

6m  11 

13:4P 

8m  29 

08:45 

10:27 

21  Jan.  59 

12s 

.7101 

15m  28  

6s  .7495 

6s  .7125 

100  gm. 

2m 

06 

21  Jan.  59 

12s  89 

12s  19 

11:00 

4m 

06 

20  gm. 

30s  83 

30s  03 

6s  0  7102 

75m 

06 

16:15 

lm  79 

lm  .7092 

12s  01 

110m 

06 

12s  .7063 

2m  70 

2m  81 

30b  00 

21  Jan.  59 

30s  69 

4m  56 

4m  73 

lm  00 

100  gm. 

In  72 

3ra  40 

3m 

2m  00 

15:31 

2m  75 

15m  30 

21  Jan.  59 

4m  .7099 

6s 

.7083 

4m  76 

30m  22 

100  gm. 

21  Jan.  59 

:  i2s 

81 

8m  78 

21  Jan.  59 

10:40 

200  go. 

30s 

80 

15m  17 

1  OH  rrm 

123  .7032 

11:05 

lm 

70 

32m  7SL 

09:17 

30s  83 

6s  .7084 

2m 

78 

6s  .7390 

lm  84 

12s  83 

4m 

77 

12s  85 

2m  84 

30s  82 

8ra 

76 

30s  74 

4m  84 

lm  81 

15n 

76 

Dote  ,Time 


Load,  Dia l 


Remarks 


— 


PLATE  42 


! 


S-* 


Minutes 


’!?£  -  Minutes 


UNIVERSITY  of  ALBERTA 

;projfct  '1 

SITE 

m 

j  DEPT  of  CIVIL  ENGINEERING 

-LAMP  LX _ Silt 

SOIL  MECHANICS  LABORATORY 

LOCATION 

CONSOLIDATIOIsjpESULTS 

HOLE 

depth 

'TECHNICIAN 

DATE 

Specific  Gravity  of  Soil  Solids  Gc  3  2.70  H*mh.  nf  smi  Saiih.  m  . 

0.613  mt 

Void  Ratio  e  (End)  3 

0.537 

Void  Ratio  e (St art)  3 

0.£08 

Void  Ratio  e  (Start  Dimensions)1 _ 

e  (End)  =  W  %lEnd)  x  Gs  Hs  =  (- 


Wt-  Soi 


ins- 


e  3  previous  e  £ 


Def'l- 


Time 

Interval 

Load  on 
Pan  (gms) 

Corr-  Dial 
Reading(ins) 

Deflection 

(ins) 

Deflection 

^s 

Void  Ratio 

e 

Press  ure 
Kg/cm1  *  T/ft* 

..  

20 

.7513 

.0434 

♦  .071 

0r  60$ 

30m 

50 

.7422 

.0343 

♦  .056 

0.593 

. 

:  .1. 

5  On 

100 

.7298 

.0219 

♦  .036 

0.573 

0.31 

15m 

200 

.7188 

.0109 

♦  .018 

.555 

0.0  0 

8m 

400 

.7069 

.0010 

/  .002 

0. 53f 

1.13 

4m 

100 

.7084 

.0005 

+  .001 

0.538 

.31 

8m 

20 

.7120 

.0041 

♦  .007 

0.544 

(.009 

50 

.7111 

.0032 

♦  .005 

0.542 

j.10 

4m 

100 

.7099 

.0020 

♦  .003 

0.540 

0.31 

9m 

200 

.7080 

.0001 

0 

0.537 

0.60 

1 47m 

400 

.7039 

.0040 

/  .007 

0.530 

l.ie 

110m 

20 

c7106 

.0027 

♦  .004 

0.T41 

15m 

100 

.7076 

.0003 

0 

0.537 

0.31 

i— - 


Compressive  Index  3 
Swelling  Pressure  3 
Pre-Cons.  Load  3 


a> 


a 

a: 


o 

> 


O-OI 


100 


pressure  Kg./cm*(Tons/ft*) 


PLATE  A  A 


' 


. 


UNIVERSITY  of  ALBERTA 

PROJECT  . 

SITE 

i  rz 

DEPT-  of  CIVIL  ENGINEERING 

SAMPLE 

SOIL  MECHANICS  LABORATORY 

CONSOLI  DAT  lONjRESULTS 

LOCATION 

HOLE 

DEPTH 

TECHNICIAN 

DATE 

Specific  Gravity  of  Soil  Solids  Gs  s 

/oid  Ratio  e  (End)  = _ 

/oid  Ratio  e(Start)  = _ 


yoid  Ratio  e  (Start  Dimensions) 


Height  of  Soil  Solids  Hs=. 


.mi 


(End)  =  W %(End)  x  Gs 


s  '  Gs  x  Area  x  £-54 ' 


ms. 


e  -  previous  e  t  — 

nc 


Time 

nt  e  rva 1 

Load  on 
Pan  (gms) 

Corr-  Dial 
Reading(ins) 

Deflection 

(ins.) 

Deflection 

Hs 

Void  Ratio 

e 

■  — ■  ■  "  ■  "g 1  » 

Press  ure 
Kg/cmJ=T/ftA 

9m 

200 

.7059 

.0020 

-  .003 ' 

0.534 

0.6Q 

15m 

400 

.7028 

.0051 

-  .008 

0.529 

1.18 

32m 

20 

.7079 

0 

0 

0^37  J 

0.069 

— 

-  ■  ■  - -t 

-1 

J 

1 

— — 

- — r 

mm — i] 

n  1  nr  - 

■  Pre-Cons.  Load  =  NIL  T/fP 


,  i  H  HI 

-I — i  i-H  Hftif  ! 


pressure  K  g/cm*(Tons/ftA) 


4-4 


■ 


173 


TABLE  XV 


CALCULATIONS  ON  CONSOLIDATION  RESULTS 


INCREMENT 
kg. /cm. 2 
(1) 

t50 

sec. 

(2) 

2H 

in. 

(3) 

ir~ 

cm.^ 

w 

2  A 

cm.  /sec . 
(5) 

2/^ 

cm.  /gm. 

(6) 

e 

(7) 

k 

cm./ sec . 
(8) 

0.069-0.16 

135 

1.01 

1.64 

0.0024 

0.00016 

0.60 

2.4  x  10"? 

0.16  -0.31 

no 

1.00 

1.61 

0.0029 

0.00013 

0.58 

2.4  x  " 

0.31  -0.60 

56 

0.99 

1.58 

0.0056 

0.000062 

0.56 

2.2  x  " 

0.60  -1.18 

3^ 

0.98 

lo4 

0.0091 

0.000034 

0.54 

2.0  x  " 

The  following  notes  apply  to  the  above  table: 

1.  The  tcQ  values  were  obtained  from  the  time  curves  which 
are  included  in  tnis  Appendix. 

2.  2H  is  the  average  thickness  of  the  sample  for  the  loading 
increment  under  consideration. 

3.  The  value  of  cv  (the  coefficient  of  consolidation),  was 
determined  using  the  equation. 


c  -  TK2 

v  "  T 

where  T  is  the  time  factor  and  is  equal  to  0.2  for  50$  consolidation 
and  the  drainage  conditions  of  the  test  performed.  (See  The  Funda¬ 
mentals  of  Soil  Mechanics  by  Donald  W.  Taylor,  fifth  printing, 
September  1951>  chapter  10 ). 

4.  The  value  of  (the  coefficient  of  compressibility)  was 
determined  using  the  equation 

av  =  ex  -  e2/p2  -  pi 


where  e  and  p  refer  to  void  ratios  and  applied  unit  loads,  respectively. 

5.  The  permeability  (k)  was  computed  using  the  equation 


k  =  cv  av  v 
1  +  e 


All  values  used  in  the  equation  are  the  average  for  the  increment 
considered. 


'  '  7.  •.  7‘  ■  •  ;  .• 


*v  1  *  7  ....  i  I  •  .  -  -  * 


f  .  .  v 


)  -•  V 


„>r  (1.  • 


174 


TABLE  XVI 
CAPILLARITY  TEST  DATA 


DATE 

TIME 

HG.  ELEV. 
Left  Right 
Cm .  Cm . 

DIFF. 

(3)-(4) 

Cm. 

DIFF. 

(5)13*6 

HEIGHT  OF 
WATER  COLUMN 
Cm. 

NEGATIVE 
PRESSURE 
Cm  of 

(1) 

(2) 

^3 

(4) 

(5) 

(6) 

(7) 

(0) 

Elevation 

of  bottom  of  specimen  was  135*5  Cm. 

Sept. 

17 

1650 

55.0 

52.7 

2.3 

31 

80 

111 

Sept. 

18 

840 

55.5 

52.2 

3*3 

45 

80 

125 

It 

930 

55.8 

52.0 

3*8 

52 

80 

132 

If 

1055 

56.2 

51.6 

4.6 

63 

79 

142 

It 

1150 

56.4 

51*3 

5*1 

69 

79 

148 

t! 

1355 

56.7 

51*1 

5*6 

76 

79 

155 

II 

1540 

57.1 

50.6 

6.5 

88 

78 

1 66 

II 

1630 

57.4 

50.3 

7*1 

97 

11 

175 

19 

1900 

57.5 

50.3 

7.2 

98 

11 

176 

II 

1000 

58.1 

49.6 

8.5 

115 

77 

192 

It 

1100 

58.5 

49.2 

9*3 

126 

it 

203 

tl 

1200 

59.0 

48.7 

10.3 

140 

11 

217 

11 

1330 

59.4 

48.3 

11.1 

151 

76 

227 

It 

i4oo 

59.8 

47.8 

12.0 

163 

75 

238 

11 

1500 

60.3 

47.4 

12.9 

176 

vi 

251 

It 

1540 

60.7 

47.O 

13*7 

187 

ti 

262 

11 

1615 

61.2 

46.6 

l4.6 

199 

74 

273 

It 

1705 

61.6 

46.2 

15*4 

210 

ll 

284 

Sept. 

20 

905 

62.0 

45.8 

16.2 

220 

11 

294 

11 

945 

62.4 

45.4 

17*0 

232 

73 

305 

11 

1110 

62.7 

45.0 

17*7 

241 

tl 

314 

11 

1310 

63.5 

44.2 

19*3 

262 

ti 

334 

It 

1600 

63.9 

43.8 

20.1 

274 

11 

346 

Sept. 

23 

900 

64.1 

43.6 

20.5 

279 

71 

350 

II 

945 

64.5 

43*2 

21.3 

290 

VI 

361 

11 

1040 

65.O 

42.8 

22.2 

302 

it 

373 

11 

1245 

65*4 

42.4 

23.0 

313 

70 

383 

11 

1540 

65*7 

42.1 

23*6 

321 

iv 

391 

11 

1630 

65*9 

41.8 

24.1 

328 

tt 

399 

tt 

2015 

66.2 

41.5 

24.7 

338 

69 

407 

Sept. 

24 

840 

69.6 

38.0 

31.6 

430 

66 

496 

II 

1505 

71*3 

36.5 

34.8 

473 

64 

537 

11 

1630 

72.4 

35*5 

36.9 

502 

63 

565 

Sept. 

25 

830 

72.8 

35*0 

37*8 

514 

62 

576 

11 

1050 

73*4 

34.5 

38.9 

530 

592 

It 

1230 

74.0 

34.0 

40.0 

544 

6l 

605 

NOTE:  Pressures  vere  applied  by  means  of  a  water-over-mercury 
U-tube.  Column  (3)  applies  to  the  water-mercury  interface  and 
column  (4)  to  the  mercury  surface  that  is  open  to  the  atmosphere. 

The  distance  from  the  water-mercury  interface  to  the  bottom  of  the 
specimen  is  given  in  column  ( 7 ) • 


TABLE  XVJ1 


DETERMINATION  OF  1*EAN  AND  MEDIAN  TE2>'JJ5RATURES 
FOR  VARIOUS  THERMOSTAT  SETTINGS 


Temp. 

°F 

Frequency 

UM2) 

Greater  Than 

Ho.  1» 

~nr 

C2; 

(3; 

(4) 

(?) 

-4.25" 

-4.75 

-5.25 

-5.75 

-6.25 

-6.75 

-7.25 

-7*75 

-8.25 

-8.75 

-9.25 

1 

4 

7 

5 

5 

8 

8 

7 

9 

8 

21 

~BT 

4 

19 

37 

29 

37 

54 

58 

52 

74 

70 

194 

-o&r 

1 

5 

12 

17 

22 

30 

38 

45 

54 

62 

83 

1.2 

6.0 

14.5 

20.5 

26.5 

36.2 
45.8 

54.2 
65.1 
74.7 

100 

,0 

Mean  Temp.  =  -7.6  F 

Median  Temp.  =  -7 • 5°F 

The  median  is  the 
quantity,  "5C#  greater 
than". 

+17.75 

+17.75 

2 

35 

2 

2.1 

17.25 

7 

121 

9 

9.4 

16.75 

3 

50 

12 

12.5 

16.25 

6 

97 

18 

18.7 

Kean  Temp.  -  +13.8°F 

15.75 

7 

108 

25 

26.0 

Median  Temp.  =  +13.8°F 

15.25 

5 

76 

30 

31.3 

14.75 

6 

86 

36 

37.5 

14.25 

3 

43 

39 

40.6 

13.75 

9 

124 

48 

50.0 

13.25 

9 

119 

57 

59.4 

12.75 

8 

102 

65 

67.7 

12.25 

12 

147 

77 

80.2 

11.75 

15 

177 

92 

95-8 

+11.25 

4 

“ST 

.JtL. 

1330 

96 

100 

+7.75 

3 

23 

3 

3.8 

7.25 

5 

36 

8 

10.0 

Mean  Temp.  =  +4.5°F 

6.75 

6 

40 

14 

17.5 

6.25 

6 

37 

20 

25.0 

The  mean  temp- 

5.75 

5 

29 

25 

31.2 

erature  is  the  total 

5.25 

4 

21 

29 

36.2 

of  column  (3)  divided 

4.75 

5 

24 

34 

42.5 

by  the  total  of 

4.25 

6 

25 

40 

50.0 

column  (2). 

3.75 

8 

30 

48 

60.0 

3.25 

10 

32 

58 

72.5 

Median  Temp.  =  +4.3°F 

2.75 

22 

“Bo" 

-52- 

356 

80 

100 

"The  temperature  given  is  that  at  the  middle  of  the  range 
covered,  i.e.  4.25  indicates  that  the  reading  fell  in  the  range 
4.0  -  4.5  degrees  Fahrenheit. 
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TABLE  X  V 1 1 1 

THERMOCOUPLE  RECORDER  CALIBRATION 
THERMOCOUPLE  TEMPERATURE 


No. 

Trial  1 

Trial  2 

(1) 

(2) 

(3) 

5 

32.5 

32.5 

6 

32.5 

32.4 

7 

32.5 

32.5 

8 

32.5 

32.4 

9 

32.5 

32.4 

10 

32.5 

32.4 

11 

32.5 

32.4 

12 

32.5 

32.4 

13 

32.5 

32.4 

15 

32.8 

32.6 

16 

32.5 

32.5 

Note :  The  calibration 
was  performed  with  the  thermo¬ 
couples  immersed  in  a  mixture 
containing  approximately  equal 
parts  of  ice  and  water. 


A  LATER  CALIBRATION 


Thermocouple 

No. 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

Ik 

15 

16 


Temperature 

33-2 

tt 

If 

It 

II 

33-3 

n 

it 

it 

n 

n 

it 

it 

33.2 

If 

It 


CALIBRATION  OF  CAPILLARY  SECTIONS 


17S 


O 

M 

8 

pa 

£ 

pa 

Q 


o 


§ 

M 

CO 

CO 


&a 


o 

t 

V 

5  to 
o  <d 

M  « 

Eh 


pa 


8 


co 


CVJ 


co 


CVJ 


ON  CO  30  ON  ON  NO  NO  LTN  ITN 

•  •  . . 

ooooooooo 


CO  On  CO  On  CO  CO  00  On  On  OncO  CO  t— 

ooooooooooooo 


ONCO  CO  ON  On  no  no  lcn  lcn  CO  OncO  OnOO  CO  CO  On  On  OncO  CO 

666666666 


H  CVJ  CO-3-  LCN-4  CO  t 


OOOOOOOOO 

•  •••••••• 

o 


qj-vqaqpw^oQooop 

H  H  H  H  CVJ  CM  CM  co_4  LCN-4  co  CVJ 


O  O  O  O  O  O  O 

H  CVJ  CO-4  LCN -4  CO 


cm  h 


ooooooooooooo 

•  •••••••*••• 

4'0  300W400000  0 
CVJ  CVJ  CVI  CO-4-  iCN-4  co  CVJ 


ITN 


ONCO  CO  CO  CO  NO  NO 

o  o  o  o  o  o  o 


ON  ONOO  ON  NO  NO  NO 

0660666 


t— 00  CO  ON  t— NO  NO 

o  o  o  o  o  o  o 


On  CO  CO  CO  CO  no  no  On  On  CO  On  no  no  no  C—  CO  CO  On  c-no  no 


OOOOOOO 

CVJ  00-4-  LCN-4-  CO  CVJ 


OOOOOOO 

OOOOOOO 

CVI  CO-4-  LCN-4-  (O  CVJ 


OOOOOOO 

CVJ  CO-4-  LTN -4-  CO  CVI 


OOOOOOO 

CVJ  CO-4-  LCN-4-  co  CJ 


OOOOOOO  OOOOOOO 

•  •••••• 

OOOOOOO 

CVI  CO-4  LCN -4  co  CVJ 


OOOOOOO 

CVJ  CO  _4  LCN -4  CO  CM 


CM 


CO 


